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VELOCITY MODIFICATION OF HI POWER SPECTRUM 

A. LAZAPJAN 1 ' 2,3 D. Pogosyan 2 

The distribution of atomic hydrogen in the Galactic plane is usually mapped using the Doppler shift of 
21 cm emission line and this causes the modification of the observed emission spectrum. We calculate the 
emission spectrum in velocity slices of data (channel maps) and derive its dependence on the statistics 
of velocity and density fields. We find that (a) if the density spectrum is steep, i.e. n < —3, the 
short-wave asymptotics of the emissivity spectrum is dominated by the velocity fluctuations; (b) the 
O velocity fluctuations make the emission spectra shallower, provided that the data slices are sufficiently 
^ \ thin. In other words, turbulent velocity creates small scale structure that can erroneously be identified 
■ as clouds. The effect of thermal velocity is very similar to the change of the effective slice thickness, but 



the difference is that while an increase of the slice thickness increases the amplitude of the signal the 
c3 \ increase of the thermal velocity leaves the measured intensities intact while washing out fluctuations. 
The contribution of fluctuations in warm HI is suppressed relative to cold component when velocity 
channels used are narrower than warm HI thermal velocity and small angular scale fluctuations are 
measured. We calculate how the spectra vary with the change of velocity slice thickness and show that 
CN ■ the observational 21 cm data is consistent with the explanation that the intensity fluctuations within 
! individual channel maps are generated by a turbulent velocity field. As the thickness of velocity slices 
; increases density fluctuations get to dominate emissivity. This allows to disentangle velocity and density 
CN ■ statistics. The application of our technique to the Galactic and SMC data reveals spectra of density and 
^ ! velocity with the power law index close to —11/3. This is a Kolmogorov index, but the explanation of 
; the spectrum appealing to the Kolmogorov-type cascade faces substantial difficulties. We generalize our 
Q\ ■ treatment for the case of a statistical study of turbulence inside individual clouds. The mathematical 
\ machinery developed is applicable to other emission lines. 

Subject headings: interstellar medium: general, structure-turbulence-radio lines; atomic 



Or 



hydrogen 

P ■ 

^ I 1. INTRODUCTION 

Galactic HI is an important component of the interstellar media (McKee & Ostriker 1977) and much 
'O \ efforts have been devoted to its studies (see Burton 1992). In particular, Crovisier & Dickey (1983) and 
5h ■ Green (1993) measured the spectrum of 21 cm intensity fluctuations in order to get a handle on the 
. 5^ , statistical properties of this media. Lazarian (1995, henceforth L95) suggested a statistical inversion 
technique for HI interferometric data (see also Lazarian 1994 a,b). Note that statistical description is a 
nearly indispensable strategy while dealing with turbulencefj. Indeterminism of turbulence implies that 
two virtually identical systems develop very different patterns of turbulent motions. The advantage of 
statistical techniques is that they extract underlying regularities of the flow and reject incidental details. 

Attempts to study interstellar turbulence date as far back as 50s (Horner 1951, Munch 1958, Wilson 
et al. 1959) and various directions of research achieved various degree of success (see reviews by Kaplan 
& Pickelner 1970, Dickman 1985, Lazarian 1992, Armstrong, Rickett & Spangler 1995). Studies of 
turbulence statistics of ionized media were successful (see Spangler & Gwinn 1990, Narayan 1992) and 
provided the information of the statistics of plasma density^ at scales 10 8 -10 15 cm. This research profited 
1 Department of Astrophysical Sciences, Princeton University, Princeton 
2 Canadian Institute for Theoretical Astrophysics, University of Toronto 
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4 Speaking about interstellar turbulence we understand it in terms of unpredictable spatio-temporal behavior of nonlin- 
ear systems (Scalo 1985, 1987) and accept that Kolmogorov picture may be too simplistic while dealing with such complex 
media (Lazarian 1999). 

incidentally the found spectrum was Kolmogorov one. 
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a lot from clear understanding of processes of scintillations and scattering achieved by theorists (see 
Goodman & Narayan 1985, Narayan & Goodman 1989). 

Deficiencies in the theoretical description have been, to our mind, the major impediments to studies of 
turbulence in the neutral part of the interstellar medium. For instance, important statistical studies of 
molecular clouds (Dickman & Kleiner 1985, Miesch & Bally 1994) has not achieved the success parallel 
to the one in scintillation studies. These problems were addressed in L95, where the quantitative study 
of deconvolution of the 3D statistics from the data averaged along the lines of sight was presented (see 
also Lazarian 1999a). 

However, the issue of velocity distortions of the observed spectrum was treated only qualitatively 
in L95, and this motivates our present study. Indeed, even in the absence of density fluctuations in 
Galactic coordinates velocity fluctuations would make the distribution of emission inhomogeneous in 
velocity data cubes. The release of the Canadian Galactic Plane survey data makes this work very 
timely. 

The present paper concentrates on the studies of HI in Galactic plane. In also addresses the issue 
of studies of individual clouds. We deliberately restrict our attention to HI, although the presented 
technique is potentially applicable to other spectral lines. The pervasive distribution of neutral hydrogen 
presents a sharp contrast to the localized distribution of molecular species and this alleviates problems 
with averaging^ (cf. Houligan Sz Scalo 1990). The remarkable constancy of the slope of the intensity 
power spectrum measured by Green (1993) and Stanimirovic et al (1998) indicates that the HI sample is 
homogeneous in a statistical sense. However, most of our results are applicable to other emission lines. 

It may seem that the Galactic plane presents simultaneously the best and the worst case for obtaining 
HI statistics. On one side, the Galactic rotation curve presents a natural distance indicator, which 
allows to roughly segregate emission from different parts of Galactic disk. On the other side, the 
random velocity is being confused with the one arising from regular rotation resulting in distance errors. 

To what extent the statistics of the underlying density can be recovered is the issue this paper deals 
with. The effects of projection and velocity distortions plague the statistical studies of turbulence in 
interstellar clouds. This paper provides a guide on how to deal with such problems. We aimed our 
paper at providing a general mathematical formalism necessary for dealing with velocity distortions 
of the observed intensity fluctuations. Modifications of the technique and its application to various 
emission/absorption lines will be done elsewhere. 

Another motivation for our study stems from recent attempts to describe the structures in the Galactic 
hydrogen in order to estimate the fluctuations of microwave polarization arising from interstellar dust. 
This contribution is extremely important in view of present-day efforts in the CMB research (see Prunet 
& Lazarian 1999, Draine & Lazarian 1999). Some of the studies, for instance, one by Sethi, Prunet and 
Bouchet (1998) attempts to relate the statistics of density observed in the velocity space and the statistics 
of magnetic field fluctuations. Our present work and a subsequent paper show that the relation between 
the two statistics is far from trivial and velocity fluctuations may create a lot of emissivity structures 
in the channel maps. 

The distortion of the density statistics in the presence of velocity fluctuations is important for un- 
derstanding of the large-scale distribution of galaxies. The corresponding literature is vast and we just 
mention a pioneering study by Kaiser (1987) and a recent review by Hamilton (1998). In cosmology the 
problem is known as the problem of 'redshift-space' corrections to the statistics of galaxy distribution 
and has been addressed predominantly under the linear regime of the gravitational instability or at the 
regime when the velocity contribution to the Fourier spectrum can be factorized, e.g. by a Maxwellian 
factor (see Hamilton (1998)). The problem that we study is much richer. We have to deal with the 
case of general, nonlinear, density fields transformed by coherent turbulent velocities. Development of 
the formalism for HI studies requires to consider a variety of spectra, including steep (e.g. Kolmogorov) 
power spectra which have never been addressed in cosmological literature. 

Our treatment of the problem is based on the description of HI statistics in terms of power spectra. 

6 A study by Stutzki et al (1998) shows that the worries about the turbulence homogeneity are probably exaggerated 
even in the case of molecular clouds. At least the wavelet analysis that they used revealed a well-defined spectrum of 
turbulence. 
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Such a description is traditional in hydrodynamics and MHD theory and it allows to relate the statistics 
of interstellar media to the statistics of idealized fluids. Moreover, spectra provide a useful intuitive 
insight to the properties of turbulence. For instance, shallow spectra correspond to more structure at 
a small scale, while steep spectra correspond to a random field dominated by large-scale fluctuations. 
Other descriptors are reviewed in Lazarian (1999b). 

An additional complication arises from variations of gas temperature. Gas thermal velocity smears 
up the turbulent velocity effects and this smearing effect depends on the temperature. In the paper we 
show that in most circumstances the intensity fluctuations are due to cold HI. 

In what follows, we formulate the problem of velocity mapping, discuss the 3D power spectrum of 
intensity that can be measured in velocity data cubes (section 2). Then we show how the 2D spectra of 
channel maps (those spectra can be directly obtained using an interferometer) depend on the thickness of 
the velocity slice (section 3). For power- law density and velocity statistics we calculate analytically and 
numerically the 2D spectra when the density is dominated by small-scale and large-scale inhomogeneities 
(section 4). We compare our predictions with the available Galactic and Small Magellanic Cloud (SMC) 
observational data in section 5, while section 6 contains a qualitative discussion of our results and 
presents an outline of the outstanding problems. We provide the summary of our results in section 7. 
To simplify our presentation we moved derivation of many results into Appendixes. They, however, have 
the value on their own. For instance, Appendix B presents asymptotics of 3D spectra in velocity space, 
while Appendix C contains derivations for correlation functions in velocity space. The modification of 
our analysis when velocity and density fields are correlated is presented in Appendix D. Appendix E 
describes a statistical treatment of emission from individual HI clouds and establishes the connection 
between this treatment and the treatment of fluctuations of intensity from HI Galactic disc that the 
main body of the paper deals with. Note, that it is this modification of the technique that is most 
important for studying turbulence in molecular clouds. Appendix F describes additional statistical 
tools, namely, ID radial spectra which can be useful measures of density and velocity. 

2. THE PROBLEM. 
2.1. The Model. 

In studying Galactic HI observers use a model Galactic rotation curve (see Kerr et al 1986) to identify 
points along the line of sight in their data cubes, in which one of the coordinates is the Doppler shifted 
velocity. The model curves reflect the averaged motion of HI in the Galactic potential or bulk deviations 
from circular motions and therefore are accurate up to turbulent velocities. 

It is easy to understand that the random velocity field makes statistics of emissivity different from 
the statistics of density. Indeed, the velocity fluctuations can make two emitting elements at different 
distances from the observer to overlap in the velocity space producing an element of doubled emissivity. 

In the present paper we address an important issue of relating the the spectrum of HI emissivity in 
velocity space xyv (v is the z-component of velocity^ to density fluctuations in real space xyz coordi- 
nates, where z is directed along the line of sight and x Sz y span the sky plane. The statistics in xyz 
coordinates is essential for understanding HI structure and processes of fragmentation (see Scalo 1985), 
while observations provide us with the statistics of velocity data cubes. 

Although the problem of relating statistics in velocity space and Galactic coordinates is quite general, 
in what follows, we limit ourselves to the case of 21 cm emission. HI observations include high precision 
measurements of Doppler shifted emission lines. For the sake of simplicity, further on we shall disregard 
self-absorption of the emitted radiation. This looks as a reasonable assumption as numerous studies 
(see Braun 1997, 1999, Higgs 1999) confirm that for many HI regions self-absorption is negligible^. 

2.2. Space- Velocity Mapping. 

7 These coordinates are also called PPV coordinates which abbreviates Position-Position- Velocity (see Vazquez- 
Semadeni 1999). 

8 In general, if absorption is localized in the form of individual blobs of dense cold material, those only alter marginally 
statistical measures on the sizes less than the size of blobs (L95). 
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The data of HI surveys is stored as data cubes with velocity v = \v z \ constituting the third, z, 
dimension. We analyze what effect the mapping from the xyz space to the xyv space has on the power 
spectrum of HI density. First, HI gas experience the regular motion induced by galactic rotation. A 
model of Galactic rotation provides a relation between the z projection (the only projection available 
through Doppler measurements) of regular velocity v r z eg and the distance along the line of sight, i.e. 
v r z eg = 0(x). If observations are made in the slices of the finite thickness the function can be expanded 
into Taylor series around a particular point (e.g. the center of data slice with v r Jf = (p(z )), v r z e9 = 
0(x o ) + 0, x (x o )(x - xo) + l/20 iXjX /(z o )(x - x )(x' - x ) + . . .. 

In addition to the regular part \ reg arising from Galactic rotation the total observed velocity of a gas 
element v o6s , also contains a random, turbulent, part u, so that v obs = V re9 + u. Therefore the mapping 
from velocity space to Galactic coordinates becomes more involved. Most interesting effects on HI 
statistics in velocity space arise from the random u part. We simplify their treatment by approximating 
the regular z-component of velocity by the linear, z only dependent term: 

z s = z s0 + A [0 )Z (z o ){z-z ) +u z ] (1) 

where z s is velocity coordinate in spectral-line data cubes, A is an arbitrary scaling coefficient that 
relates velocity and data cube length units. Naturally, a linear mapping given by eq. ([!]) is only correct 
if the nonlinear terms in the expansion of are small ( in other words, we consider the region of space 
where coherent flow can be described as flow of a constant shear) and the scales of interest in image 
plane are sufficiently short that x,y dependence of the shear can be ignored []. 

To describe turbulence in Galactic disc we find advantageous to use parameter f(zo) = 0^(2o) _1 and 
put A = f. Using vector notations we describe the effect of a random velocity u as a linear map from 
Galactic frame x to the velocity-space coordinates x s 

x s = x - /(x ) (u • z) z , (2) 

Here z is a unit vector along the line of sight, x is the actual Galactic coordinate of the emitting HI 
element and x s is the position calculated using the Galactic rotation curve (see Blitz & Spergel 1991). 
We shall note that the values of / vary from slice to slice in accordance with the Galactic rotation 
curve. In our Galaxy / is of order of 100 pc/ (km/s), although it depends on the direction and depth of 
observations. 

A limiting case of the mapping corresponds to the study of individual HI clouds. We discuss this case 
in Appendix E and present a generalized form of the map applicable for the zero shear, i.e. f^ 1 — > 0. 

2.3. Spectrum in the Velocity Space 

The relevant expression for the spectrum of density in velocity space has been derived by Scoccimarro, 
Couchman & Frieman (1999). Let us write the Fourier component of the density p(x s ) in velocity space 
as the sum over positions x" of all hydrogen atoms 

p s (k) = £ e* kx * = J2 e^ a - f( - k ^ ua ^ , (3) 

atoms atoms 

In a coarse-grained fluid description of HI atom density in real space, the summation over all atoms 
is replaced by a continuous integral in the xyz (real space) coordinates: J2 a toms =^ I ^ 3 xp(x) with gas 
density p(x). Velocities of individual atoms are averaged^] over the thermal distribution (exp[— if(k ■ 

9 It is obvious that if the velocity of Galactic rotation is V max then <j>, Z z(zo)(z — Zo)/<f>,z(zo) ~ (SV/V max ), where SV is 
equal or less than the thickness of a velocity slice. For sufficiently small slice thickness the factor (5V/V max ) is small and 
the contribution of the nonlinear term is negligible. It will be clear from the discussion below that studying turbulence 
at a particular scale one should account for the non-linearity of <f> on this scale. Similarly, fluctuations are observed over 
a limited area of SxSy plane and one can use <p(zo, xq, yo) to characterize the regular velocity if variations of <fi on the 
distances studied are small. This variation decreases with a distance to the slice. It is easy to see that this requirements 
for the map to be linear are easy to fulfill. 

10 The averaging is similar to that presented by eq. (]To|) . 
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z)(u a -z)]) thermal = exp [— ifk z u z (x)] exp [— f 2 k 2 v^/2] , where u(x) is now the peculiar turbulent velocity of 

/ r X 1 / 2 

the fluid element at position x. Thermal velocity with one dimensional mean-square vt = 0.75 ( yo^k ) 
km/s is incoherent and acts as a smoothing along the velocity coordinate, limiting the resolution in 
velocity space. Thus we obtain for density Fourier amplitude 

p s (k) = e - /2fc ^ /2 y d 3 xp(x)e- i/fc ^ z(x) e ik ' x . (4) 
The variance of Fourier amplitudes in velocity space is 

(p s (k)^(k')> = e -f 2 (kl+K 2 K/2 [ rf s x j d 3 x / e acx-k'.x' ( p ( x ) p ( x /) e -i/[W*)-fc>*(*')]) . (5) 

Henceforth we use angular brackets (...) to denote ensemble averaging over turbulence realizations. We 
assume statistical homogeneity of the ensemble in the Galactic frame, i.e. we assume that the average 
of any expression with two spatial coordinates depends only on the vector separation between these two 
points. The density Fourier modes in velocity space are then uncorrelated and can be described by the 
power spectrum^ P s (k) 

(p s (kK(k')> = P,(k)5(k - V)e-^ ( 6 ) 
P s (k) = yd 3 re ik - r H(k,r), r = x - x/ , (7) 

where the kernel is 

S(k,r) = (e^ (Mz(x) - Uz(x/)) p(x)p(x/)) . (8) 

Looking at eq. @ one easily notices that although both velocity and density contribute to P s (k) their 
functional dependences are different. This opens a perspective of separating the velocity and density 
contributions to the emissivity statistics that we explore in sections 4 and 5. 

Eq. (||) assumes a plane-parallel approximation when one does not distinguish between the radial 
nature of the line of sight and the Cartesian ^-direction. The formula is accurate when the HI region is 
observed from a distance that is large compared to its thickness C s i. The statistical analysis when the 
converging geometry of the lines of sight is accounted for is more involved (see Lazarian 1994b) and will 
be discussed elsewhere. 

In what follows we limit ourselves with the discussion of basics of the space- velocity mapping and the 
role of coherent turbulent velocity. It is easy to see that the thermal velocity effect in z s direction is 
similar to the effect of the finite velocity resolution of the telescope in the x and y directions (see L95). 
To make our formulae more comprehensible we defined the spectrum P s (k) to contain only correlated 
turbulent contribution while keeping thermal velocity factor separately. We can neglect thermal effect 
altogether if we study supersonic turbulence in velocity slices 5V > vt- We, however, will remember 
about the high k z truncation of our spectrum by thermal motions when necessary (see Chapter 3). 

To proceed at this stage we make the following assumptions about the properties of the random fields 
p(x) and u(x) . First of all, we assume that the turbulent velocity is uncorrelated with the density in 
Galactic coordinates, and therefore (e l ^"p(x)p(x/)) = (e l -^")£(r), where the density correlation function 
is assumed to be isotropic in xyz-sp&ce 

£(r)=£(r) = (p(x)p(x + r)> . (9) 

We discuss the corrections due to velocity-density correlations in Appendix D and show that they 
are of limited importance for Lognormal model of density fluctuations. Future research will test our 
assumption for the distributions produced by numerical simulations. Interstellar velocity and density 
fields can be anisotropic and we discuss this in sections 6.2 and 6.3. 

11 Here we deal with small scale coherent structures and do not deal with the large scale ones, e.g. spiral arms. 
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We do not have to make any further assumptions about density statistics. The turbulent velocity we 
assume to be Gaussian, then 



A«; = W;(x) - Uj(x/) . (10) 



To simplify our analysis we also restrict ourselves to the case of locally isotropic turbulent velocity 
fields^. Thus the structure tensor (AuiAuj), which fully describes an isotropic vector field can be 
expressed via longitudinal Dll and transverse D^n components (Monin & Yaglom 1972) 

(A Ul A Uj ) = (D LL (r)-D NN (r))^ + D NN (r)5 lk , (11) 

where 8ik equals 1 for i = k and zero otherwise. We define z-projection of the velocity structure function 

as 

D z (r) = (AuiAuj)ziZj = D NN (r) + [D LL (r) - D NN (r)] cos 2 9 , cosd = r ■ z (12) 
Substituting eqs. ( jTT|) and ( jT0"D into eq.(^) we obtain: 

pe z D z (ry 



P s (k) = J tfW k - r £(r)exp 



(13) 



It is easy to see that the power spectrum above is anisotropic, with velocity mapping having most effect 
on the modes parallel to the line of sight (k z = k) and not affecting the perpendicular ones (k z = 0). 

2.4. Velocity field 

If the velocity field is solenoidalQ then according to Monin & Yaglom (1972) Dll and Dnn are related 
via a simple relation: 

D NN (r)=D L L(r) + r -^-D L L(r) . (14) 
2 ar 

Observational data that we discuss in section 5 corresponds to a power law. This is suggestive that 
the underlying statistics is also a power-law. For the power-law statistics velocity eq. ([TJ]) provides 
D LL = 2/(m + 2)D NN = Cr m . 
For potential fields Monin & Yaglom (1972) show that 

D L L(r) = D NN (r) + r -^-D NN (r) , (15) 
2 ar 

thus for power-law structure functions Dll = (1 + m/2)Djs[N and again Dll ~ r ™ ''■ 
In general, the velocity field will have both potential and solenoidal components and this introduces 

an uncertainty of the order unity in the coefficients relating Dll and Dnn- This does not change much 

in our results, and therefore, for the sake of simplicity, doing calculations we shall assume that the 

velocity field is solenoidal, which probably is not too far from the truth. 
The convergence of integrals restricts m to the range < m < 2. The value m = 2/3 is distinct as it 

corresponds to the Kolmogorov turbulence, but in Appendix B we consider also cases m = 1/2, 1. The 

power spectrum in velocity space is then 



P s (k) = J tfW k - r £(r)exp 



(fc 2 A) 2 ^(r/A) l 
2 

J D z (r/A) = (r/A) m (l + m/2(l-cos 2 ^)) , (17) 

12 The issue of density anisotropy is analyzed in section 6.2 and the data reduction technique dealing with anisotropies 
is discussed in L95. 

13 Solenoidality of the velocity field follows for the incompressible fluid from the continuity equation. 



7 



where the turbulence scale in velocity space is 



2f~1 2 — m 



A = [f'C 

This scale depends on the Galactic rotation curve, direction of observations and turbulence intensity. 
The scale A corresponds to the scale over which random velocities map density fluctuations. Physically 
this is the scale at which the velocity dispersion ~ CX m becomes equal to the squared difference of the 
regular velocities determined by Galactic rotation (i.e. f~ 2 X 2 ). 

While we discuss the properties of the 3D spectrum in velocity space in Appendix B, the 2D spectrum 
of HI intensity fluctuations available via observations is discussed in the next section. 

3. THICK AND THIN SLICING OF DATA CUBES 

It is easy to see that a radiointerferometer gets the 2D spectrum of 21 cm intensity (see L95). Therefore 
it is very important to relate this 2D spectrum to the underlying 3D statistics of turbulence. The 
velocity resolution of an individual radiointerferometer channel represents the thinnest slicing of the 
data available. In practice, several channels can be combined to provide a thicker slice (see Green 
1993). 

Henceforth we use capital letters R and K to denote 2D quantities and reserve small letters r and 
k for 3D quantities. In particular, R is the vector in two dimensional plane orthogonal to the line of 
sight, and the distance in xyv space is r 2 = ~R 2 + z 2 , where z s is the z-coordinate in velocity space. 

An important property of 21 cm emission is that it does not depend on gas temperature (see Spitzer 
1978). For negligible self-absorption (see discussion in section 2.1) the HI intensity at the point R, 
observed in velocity slice 5V, is proportional to the integral of the density over the slice thickness 



KB) 



J dz' sP (R, z' s ) W e 



(19) 



The width of the normalized (/ dx s W e = 1) experimental window function W e , £ s \, represents the 
velocity slice thickness; z s is the value of the slice central velocity. We expressed both quantities 
in length units, i.e. £ s \ ~ fSV, keeping in mind the coordinate transformation (||]). The shape of 
the window function depends on the sensitivity profile and the width of the individual interferometer 
channel and on how the channels are combined. In idealized case of uniform sensitivity one can consider 
step-like window W e (Az/£ s i) = l/£ s \ if \Az\ < £ s \/2 and W e = otherwise. Intensity 1(B) implicitly 
carries dependence on z s , which we omit for brevity. 
The normalized correlation function of intensity fluctuations is 



(<5/(R 1 )5/(R 2 )) 



tfe^dz^ (<fp(Ri, z' sl )5p(B 2 , z' s2 ))W t 
dz' s2 dz' sl £ s (Ri - R 2 , z' sl - z' s2 ) W e 



Z s i — z\ 



si 



£ si 



We 



Zs2 — Z s2 

£ si 



Z s\ Z sl 



We 



— Z s2 

£ s i 



(20) 



where 51 and Sp mean, respectively, variations of intensity and density. The density correlation function 
in the velocity space £ s is the Fourier transform of the turbulence spectrum P s (K,k z ) (see eq. (|T3|)) 
smoothed by the thermal factor 



(2tt 



-3 



d 3 ke- 



ik-r. 



P s (k)e 



(21) 



The properties of the correlation function are discussed in Appendix |C[ 

The sum of the squared imaginary and real parts of the interferometer visibility function is proportional 
to the two dimensional spectrum of intensity fluctuations (see L95), which is 



P 2 (K) = J rf 2 Re iKR (5/(R 1 )(5/(R 1 + R)) 



(22) 
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Substituting eqs. (20) and fl2~l]) into eq. (|22"|) after simple calculations similar to those in Appendix E 



we derive the expression for -P2(K) 

P 2 (K)| £ ~— / dk z P s (K,k z )W 2 (k z C sl )e^ k ^ , (23) 

where \c reflects the dependence on the slice thickness. W e (k z C s \) is the Fourier transform of the 
experimental window function. For instance W 2 (k z C s \) = 2 [1 — cos(k z £ s \)} /{k z C s \} 2 for the step-like 
window. 

The effective filter W{k z C s \, k z vr) acting on 3D turbulence spectrum in eq. (|23"D is the product of 
experimental window function and the thermal term 



W(k z C sh k z v T ) = W'^k z C sl )e-^T . (24) 

This demonstrates that thermal dispersion can be treated as part of the experimental velocity window. 
The effect of velocity slicing 5v is similar to the effect of vt in the sense that they both limit the k z 
resolution. Increase of smears the fluctuations. Thus we may expect the fluctuations from colder, 
vt < Sv, parts of the medium to dominate over the contribution from warmer, vt > Sv, parts. The 
exact criterion when this is true is given below. 

At small k z we have W oc 1 — k 2 (Wg"(0)C 2 i + 2f 2 v T )/2, thus the characteristic scale of the window 
can be defined by an effective slice thickness 

C = 2(W 2 "(0)£ 2 l + 2f 2 v 2 T ) 1/2 = 2f(W 2 "(0)5v 2 + 2v 2 ) l/2 . (25) 

This is the customary definition of a characteristic scale as twice the value of k^ 1 at which the filter 
value drops by a factor of two of its maximum. This approximately corresponds to the full width at 
half maximum of the window in velocity space and is the effective width of the velocity slice. In what 
follows we denote the window function given by eq. ( |24"D by W(k z C) and disregard its dependence on vt 
when it is possible. In section 4.3, we, however, will have to treat the dependence of C on Vt in order 
to determine the relative contribution of the warm and cold HI. 

For step-like experimental window, W 2 "(0) = 1/6. In this case thermal velocity dispersion has much 
more profound effect than the averaging over the velocity channel of similar thickness. This is due to 
relatively weak cutoff in Fourier space provided by step-like averaging in velocity slice. 

As long as the exact shape of the high k z cutoff by the window function is not important, using 
effective slice thickness fully accounts for the effect of finite gas temperature. From now on we shall 
mean 'effective' everywhere we speak about slice thickness. 

The thinner is the velocity slice C, the closer to unity is the window W over larger range of k z and 
more 3D modes contribute to 2D spectrum. 

We shall call a slice thin, when the W(k z C) can be approximated by unity over all k z of interest and 
the one of thick, when it is important that at high k z the function W(k z C) decreases . Therefore when 
the slice is thin eq. ( P3"D givesQ 



P 2 (K)\ t ^±-Jdk z P s (K,k z ) , (26) 

where \ t denotes that the slice is thin while in the limit of a thick slice the exact expression for the 
window function W(kC) (i.e. eq. (0)) should be used. Obviously enough, whether the slice is thin or 
thick depends not only on £ but also on the properties of P S (K, k z ). The exact criteria for slices to be 
thin and thick are discussed in the next section. 

14 Thc intuitive notion of the thin slice is hardly associated with the integral over spectrum. Note, however, that the 
integral is taken in fc-space and, for instance, a thin sheet in a real space corresponds to an integral in fc-spacc. A different 
way of understanding eq. (|2^) is to consider randomly moving 3D waves that produce a pattern corresponding to the 
wavevector K on their intersection of the plane. Obviously enough, various 3D waves with k = (K, k z ) contribute to the 
pattern. Integration over k z in eq. (|26|) reflects this. 
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One expects different behavior of resulting 2D spectrum P2 depending on a slice being thin or thick. To 
clarify this distinction, let us consider the density power-law 3D spectrum P S (K, k z ) oc (|K| 2 +fc 2 : ) n//2 , n < 
— 1 that is not subject to a velocity mapping. For this simple case -P S (K, k z ) is nearly constant if k z < |K| 
and oc k z if k z > |K|. This fall-off at high fcj (n is negative) amounts to a cut-off of the integral 



(K,fc z ) 



-l K l 



id "": 



K,k z 



at fc z ~ |K| (assuming, as we did, n < —1). Correspondingly, 



dk z P s 

if £|K| <C 1 the window function is unimportant and the slice is thin, while the slice is thick when 
£|K| ^> 1. In other words, the transition between thick and thin slices occurs when we consider 
wavelengths longer or shorter than slice thickness. For the thick slice, the contribution to the integral 
^3|) for the case considered mostly comes from the k z ~ modes i.e. 



P 2 (K)\ T ^ P s (K,k z ) J dk z W(k z C) /2tt ^ C-'P^O) ,0<k z <l/£, 



(27) 



where \j- denotes that the slice is thick, and this result corresponds to the finding in L95. The result 
can be easily understood. For sufficiently steep spectra (n < —1) integration over z corresponds to 
choosing k z = mode in the Fourier space. Thus the 2D intensity spectrum and the underlying 3D 
density spectrum with k z = should be identical. 

Thus, the asymptotics of the 2D power spectrum produced by integrating density fluctuations along 
a line of sight are 

f IKI" £IKI :» 1 
ft(K) oc I j^V, *|K|>1 (28) 

However, the rule of transition between thin and thick slices gets modified when we consider the 
emissivity power spectrum in velocity space. In the presence of velocity mapping (see eq. (|2|)) an 
additional scaled] A (see eq. (|7])) appears. This alters the wavenumber |K| at which the transition 
between the thin and thick slices occurs. 



The thin-slice expression for 2D spectrum in the velocity space follows from eq. (|26| ) with 3D kernel 
given by eq. (p~6| ) 



P*(K) 



(2vr) 1 /2A 



d 3 r e 



iKR 



£(r) 



£*(r/A)i/2 



exp 



2A 2 ^(r/A) 



(29) 



The exact relation between -P2(K) and the underlying velocity and density statistics depends on 
whether the density spectrum of HI is shallow or steep. In the next section we discuss these two cases 
separately. 

4. TWO DIMENSIONAL SPECTRA 

Below we consider the distribution of density dominated (a) by fluctuations at small scales and (b) by 
fluctuations at large scales. As we mentioned earlier, we limit our discussion to the power-law spectra 
both because observations suggest power-law dependences (see Green 1993, Stanimirocic et al 1999). If 
future research suggests a different dependences for interstellar statistics, our formulae can be used to 
fit the data. We show that in terms of power-law spectra case (a) corresponds to a power-law index 
n > — 3 and case (b) to n < —3. 

4.1. Short-wave dominated spectrum of density field 

To describe the statistical properties of the density fluctuations dominated by short wavelengths, we 
use power-law correlation functions of over- density: 



£(r) = <Pr[l+[f) J. 7>0 . (30) 

5 If the underlying 3D spectrum has several built-in scales, detailed analysis can distinguish between several intermediate 
asymptotics. 
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The power-law part of the correlation corresponds to the 3D power-law spectrum 

P(k)ock n , n = 7-3>-3 . (31) 

The unity term in eq. (|30| ) reflects the fact that density has non-zero mean value. After Fourier transfor- 
mation this term leads to a delta function term at k = added to the power spectrum. The correlation 
scale ro gives us the second parameter of the problem (the first one being A). 

Two terms in the expression for correlation function lead to the correspondent split of the power 
spectrum in velocity space into two parts. We show in the Appendix B that the 3D power spectrum in 
velocity space can be presented as 

P s (\K\,k z ) = (p) 2 [P V (\K\, k z ) + P P (\K\,k z )} , (32) 

where we introduced a notation 

P v (\K\,k z )= J d 3 re ikr exp 

P,(|K|,A*) = jd 3 re ikr (r /r) 7 exp 

In eq. ( |32|) P v (k) is the emissivity spectrum in velocity space arising from random velocity fluctuations 
in an incompressible fluid and in terms of mathematics arises from the unity term (mean density) in 
eq. fl3~0f ), while P p (k) describes the spectrum of emissivity fluctuations arising from both velocity and 
density fluctuations. Formally -P p (k) arises from velocity mapping acting upon the density fluctuations 
given by (r/ro) 7 term in eq. flSTf ). As expected, the P v term has a ^-function behavior at k = 0, but in 
velocity space it is also non- vanishing for k > 0. As we mentioned earlier even uniformly distributed in 
space, but turbulently moving incompressible fluid produce fluctuating emission (given by P v (k)) when 
observed in velocity slice of finite thickness. 

We derive the asymptotics for P v and P p in Appendix B. The main result is that the velocity mapping 
leads to the rapid cutoff of modes with large k z component of the wave vector. Both P p and P v falls 
sharply (cx k~ 2a l m where a is 3 for P v or — n for P p ) for k z ^> X^ 1 (kX) m ^ 2 , k 2 = K 2 + A; 2 . This implicitly 
defines the cutoff scale k z which can be approximated as k z = A _1 (l + |K|A) m//2 . 

We will see that P v and P p emerge, although in another combination, when the density spectrum is 
long-wave dominated. Thus, below, we formulate a general criterion for the slice to be thin or thick in 
the presence of velocity mapping. 

As we discuss in section 5, for scales of observational interests |K|A 3> 1. For these modes k z ~ 
A _1 (|K|A) m / 2 < |K|, since m < 2. This strongly affects the integration in eq. (p3|) . Indeed, one can 
assume that the window function W(Ck z ) ~ 1 as long as Ck c z < 1. In particular, the slice remains thin 
for much larger |K| than one would expect from the example of power-law spectrum unmodified by the 
velocity mapping (see eq. (|28|)). Now, given |K|A ^> 1, the transition between the thin to thick slice 
regimes is 

thin slice : C/X ■ (|K|A) m/2 < 1 

thick slice : C/X ■ (|K|A) m/2 > 1 (35) 

The thin-slice asymptotics is achieved when the first condition is satisfied due to very steep cutoff of 3D 
spectrum at high k z . In other words, due to a rapid decrease of 3D spectrum P s the window function can 
be ignored over the whole range of k z . The last inequality essentially gives us the necessary condition for 
the variations of the window function to be important. To appreciate the effect of the velocity mapping 
let us rewrite it as £|K| > (|K|A) 1_m / 2 ^> 1. This is to be compared with thick slice condition £|K| > 1 
which would take effect if there is no velocity term, A — > (combined rule £|K| > max(l, (\K\X) l ~ m l 2 ) 
is valid for any A). For the wavenumber |K| = (10 pc)" 1 , typical velocity turbulence scale A = 3.5 kpc 



(k z X) 2 D z (r/X) 



(33) 



(k z X) 2 D z (r/X) 



(34) 
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and 



(see Section 5) and m = 2/3 we thus need C > 500 pc for slice to be thick, which is fifty times larger 
than Eq. (|28|) suggests. Indeed, according to Eq. (|2~8"D C > |K| _1 = 10 pc is sufficient for the slice to be 
thick if the velocity field is absent. 

The separation between thin and thick regimes is intuitively clear if rewritten in physical units. Sub- 
stituting A from eq. ([18]) and expressing slice thickness through the width of window function, which 
for v t = is proportional to the width of the interferometer channel SV, £ = fSV we obtain 

thin slice: C\K\- m > 5V 2 

thick slice: C\K\- m <^5V 2 (36) 

In other words, if the velocity dispersion Cr m on the scale |K| _1 is larger than the squared width of the 
channel^ (in velocity units) the slice is thin. If the opposite is true - the slice is thick. We stress, that 
the critical wavelength depends only on the width of the interferometer channel and the amplitude of 
velocity turbulence, but not on the direction of observations or the actual physical thickness of the slice 
involved (al long as |K|A ^> 1). To observe the transitions between the two regimes one can either vary 
5V or |K|. We show in Appendix E that the criterion ( j3lf) is quite general and applicable not only to 
the Galaxy but to individual clouds and external galaxies. 

The 2D spectrum of intensity fluctuations in a velocity slice of data can be obtained using eq. (p3[). 
It follows from eq. (|32|) that 

P 2 (|K|) = (p) 2 [P 2 ,(|K|) + P 2p (|K|)] , (37) 

where 

P 2v (K)\ c ~ f°° dk z P v (K,k z )W(k z C) , (38) 

Z7T J ~oo 

P 2p (K) \c ~^- [°° dk z P p (K, k z ) W(k z C) , (39) 

The calculation of appropriate integrals is straightforward using the asymptotics for 3D spectra P v and 
P p obtained in Appendix B. In Table 1 we present asymptotical formulae for P 2p and P 2v spectra in a 
thin and a thick slice regimes under the assumption |K|A ^> 1. P 2p and P 2v also appear in the case of 
long-wave dominated density fluctuations and therefore our results are quite general. 

In the regime of a thick slice Table p] shows that the power-law indexes of the 2D spectrum and the 
underlying 3D density spectrum coincide. This is the regime considered in L95. If the slice is thin, the 
2D spectrum index becomes shallower, i.e. n + m/2. However this spectrum is still steeper than the 
spectrum of density in a thin slice in the absence of velocity modulations. Indeed, in the latter case 
eq. ( ^8|) predicts the index n + 1 > n + m/2 for m < 2. This thin- slice behavior can be immediately 
found from the integral representation given by eq. (P9"D- As can be easily shown, at high |K| the 
exponential factor can be set to unity, so 

p ^i^S^m^ (40) 

An extra factor D l J 2 (r) oc C x l 2 r m l 2 is responsible for the change of slope of the spectrum due to velocity 
effect. 

Note that in a thin slice P 2v behaves as a density part P 2p taken with n = —3. Comparing expression 
for P 2p rewritten as P 2p = 5 , nm (|K|r ) 3+ri (|K|A) _3+m '/ 2 with the formula for P 2v shows that for short- 
wave dominated spectra n > —3 the density part is always dominant over P 2v at wavelength shorter 

16 Here, for the sake of simplicity, we talk about measurements for which the thickness of the slice is determined by the 
velocity channel of the telescope. In reality, data from many channels can be added to produce a thicker slice (see section 
6). Moreover, in most cases vt contribution is important in eq. ( |25| ) and in more precisely SV = 0.5(W / e 2 "(0)<5w 2 + 2v^) 1 / 2 . 
Generalizations of our criterion to these cases are trivial. 
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thick slice 


thin slice 




C/X- (|K|A) m / 2 > 1 ; 


C/X- (|K|A) m / 2 < 1 


A- 2 P 2p (K) : 


A„-(A//:)-(ro/A) 3+ "-(|K|Ar; 


5„ m (r /A) 3+ "-(|K|A)" +m / 2 


A- 2 P 2 ,(K) : 


P m -(A/£) 3 -(|K|A)- 3 -™; 


5_ 3m -(|K|A)- 3+m / 2 



Table 1 

ASYMPTOTICS OF THE 2D SPECTRUM IN THE thin AND thick VELOCITY SLICES FOR |K|A » 1. NUMERICAL 

CONSTANTS A n , B m AND S nm ARE GIVEN IN APPENDIX A. IN APPENDIX E WE SHOW THAT THESE 
RESULTS ARE APPLICABLE TO INDIVIDUAL CLOUDS PROVIDED THAT THE CLOUD SIZE IS USED INSTEAD OF 

A. 



than the density correlation length |K|r > 1. The density term P 2p is also dominant for thick slices. 
Indeed, 

^(|K|r ) 3+ " (£f(|K|Ar »1 , (41) 

as shallow density with n > —3 and the terms in the second bracket ^> 1 for a thick slice. In this regime 
the velocity part has an extra (A/£) 2 factor and extra damping (|K|A)~ m . 

The magnitude of P 2v (K, k z ) is more tangible than one might suggest in an analogy with eq. (p7|). The 
reason is that P 2v (K, k z ) is a growing function of k z up to k' z = A _1 max 1, (|K|A) m / 2 (see Table 2 in 

Appendix B). This is the case where the exact shape of the fall-off of the window W(Ck z ) is important. 
Typically if W(Ck z ) falls fast at k z > C^ 1 we have 

A- 2 P 2 „(|K|) cx \£'% M V^ oc (A/£) 3 (|K|A)- 3 - m (42) 

Peculiar is the case of step-like uniform window, which falls off only as l/(k z C) 2 , then 

A- 2 P 2 „(|K|) cx / o (|K|A)W2 d(Xk z ) { [^Z n {k ' c)2 oc (A/£) 2 (|K|A)- 3 -/ 2 (43) 

Our numerical results have been obtained for this case. 

Details of the transition from thin to thick slicing as £ increases require numerical study. The results, 
plotted in Figure 1, confirm the accuracy of our analytical estimates. Our calculations show how with 
the increase of slice thickness the spectral index of P 2p steepens from the value n + m/2, characteristic 
to thin slice, towards the value n measured in thick slices. The transition scale between these two 
regimes is significantly increased as the result of velocity fluctuations. For thin slices and large K the 
P 2p spectrum depends on both velocity and density fluctuations. The pure velocity term P 2v is always 
subdominant at short wavelengths. 

All in all, on short wavelengths |K|A ^> 1 thick velocity slicing of 21 cm data cubes provides the 
spectrum |K| n , while thin slicing provides |K| n+m//2 . Therefore if density statistics is dominated by 
small scale fluctuations, interferometers can provide us both with density and velocity statistics; the 
density spectral index being obtained from thick slices, and combining this with the thin slice data one 
can obtain the velocity spectral index m. 
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Fig. 1. — Two dimensional power spectrum measured in a slice of variable thickness L — C/X (upper curves correspond 
to smaller Li). Left panel shows f|j(|K|) in the absence of velocity mapping. The middle and the right panels correspond 
to densit y Pp _ n and velocity Pi v terms. The upper dashed lines show the thin slice power-law asymptotics described 
in Table |C3| ; for comparison the short-wave Pi p slope is also shown in the left panel. The curve of enhanced weight 
corresponds to L = 1. The thick slice asymptotics are accurate below the dotted line. 



4.2. Long-wave dominated spectrum of density field 

Our analysis above does not cover 7 < and n < —3. This is the regime dominated by large-scale 
density fluctuations that we consider in this section. 

Kolmogorov turbulence is often assumed to be valid not only for the velocity but also for the density^]. 
This type of turbulence corresponds to 7 = —2/3 = — m, n = —11/3. For 7 < we formally substitute 
in Eq. (pD 

£(r)=d(oo)[l-d(r)/d(oo)] , (44) 
where d(r) is a structure function^ of density 

rf(r) = ((p(x + r)-p(x)) 2 ) . (45) 

Real world structure functions do not grow infinitely and therefore we have to introduce a cutoff at 
some large scale. If the cut-off happens at r c then for the long-wave dominated turbulence (7 < 0) 

d(r) = d(oo) r _ 7 , (46) 
r t + r c 



^) = ^-Z^ • (^) 



and the correlation function is 

d(oo) 

2 r~ 7 + rc 

For simplicity^, further on we assume that the structure functions of velocity and density have the 
cut-off at the same value of r c . Below we show that our final results marginally depend on a particular 

17 This is exactly true for the density of a passive contaminant transported by the Kolmogorov velocity field. It is not 
a priori clear that such a treatment is applicable to HI. 

18 We could start from the structure functions from the very beginning but we preferred to use correlation functions to 
make our presentation more uniform. 

19 This is a consistent assumption for the Kolmogorov density turbulence. 
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value of r c as long as r c is much larger than the turbulence scales under study, i.e. |K|r c ^> 1. Similarly, 
it is possible to show that in the presence of two different cut-off scales r ic < r 2c the results do not 
depend on either scale provided that |K|ri c ^> 1. 

Observations indicate that the amplitude of random velocity is high enough, so that A > r c (see 
section 5). Thus we are predominantly interested in the short-wave regime |K| _1 < r c < X. 

For sufficiently small r Cr c eq. ( |4"T| ) gives 

e(r)«^(oo)(l-[r/rD (48) 



Note that in eq. (|4q) r c plays the role of Tq in eq. (p0|). Both r c and r are the scales at which most of 
the energy of the field is stored. It is straightforward to show that at short waves |K|r c > 1 eq. fl4*B|) 
entails 

P 2 (\K\)^^d(oo)[P 2v + P 2p ] , (49) 

where P 2p coincides with the definition in section 4.1 if a change r c — > is made. P 2v arises from the 
unity term in the definition of correlation functions and therefore is identical to that in section 4.1. 



Eq. ([49]) is an analog of eq. ( p7\) for the case of a density field dominated by large-scale fluctuations. 
The comparison between the two reveals a substantial difference between the two regimes. Note, that 
by construction P 2p is the only part which carries the information about the index of underlying density 
statistics^} The other term, namely, P 2v arises purely from velocity. Notice, that in this case the 
variation of the amplitude d(oo) of the density fluctuations in xyz space affects P 2v and P 2p (see eq. (jiTf)) 
in the same way. Asymptotics of P 2v and P 2p are given by Table 1 provided that r c is used instead of 
ro in the expression for P 2p . 

The feasibility of extracting the density statistics from the observable emissivity spectrum depends 
on the thickness £ of the data slice, which governs the relative magnitude of P 2v and P 2p terms. In a 
thin slice, i.e when (C/X) ■ (|K|A) m//2 > 1, a comparison between P 2v and P 2p is similar to that carried 
out in the previous section. It shows that for n < —3 P 2v always dominates if |K|r c > 1. In thick slice, 



eq.([12|, P 2v has the same spectral index as P 2p but its amplitude differs by (X/C) 2 . P 2p dominates when 

C/X > (ro/A)™/ 2 (50) 

This condition does not qualitatively change when thermal line broadening is taken into account. In 
general, for C ^> A the density spectrum should be observable. 

It is clear that sufficiently weak velocity turbulence cannot distort the density statistics. Our formulae 
confirm this intuitive prediction. As the level of velocity turbulence decreases so does the scale A given 
by eq. ([i~8|). As the result C for very weak turbulence becomes much greater than A and therefore the 
density statistics reveals itself. 

Thus we can summarize that when C < X the velocity term dominates, producing universal spectral 
slope independent of the underlying density slope n. This means that for the thick slices the expected 
slope is —3 — m/2 and for thin slices the slope is —3 + m/2. If the turbulent velocity as index m = 2/3 ( 
the Kolmogorov index), the slope is equal to —8/3 if the slice is thin and steepens to —10/3 for thicker 
slices. 

The effect of velocity mapping on the density power spectrum is illustrated in Figure |2|. An additional 
asymptotics, namely K~ 4 ^ 3 , arises from the evident expansion of eq. ( |4"TD for r 3> r c . In this regime 
density fluctuations dominate those of intensity. 

There may be an element of confusion since "Kolmogorov" density fluctuations in the absence of ve- 
locity mapping will also produce n = —8/3 spectrum in a thin slice. However, as Figure |2| demonstrates, 
in the regime when random velocities are large, it is the velocity term that dominates for large |K|. The 

20 The amplitude of P2 P for n < —3 is negative and this compensates for a minus sign in eq.(^). 

21 We emphasize that the name does not imply the actual connection to the picture of the Kolmogorov energy cascade, 
but only to the particular power index. 
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Fig. 2. — Two dimensional power spectrum for "Kolmogorov" turbulence (which does not mean more that the 
statistics of density and velocity have the Kolmogorov indexes m = —7 = 2/3) measured in a slice of variable thickness 
L = £/ \ (upper curves correspond to smaller Li). The ratio X/r c = 50 is fixed. Upper left panel shows P2GKQ when the 
velocity mapping is absent. All quantities are expressed in units of r c for this panel. The bottom left panel presents the 
density P 2v term arising from velocity fluctuations only. The upper right panel shows P 2v + P 2p when both velocity and 
density are "Kolmogorov". The dashed lines show power-law asymptotics described in detail in Section 4. Heavy line 
corresponds to the slice thickness C = 0.1 A. The bottom right panel repeats the previous calculations but for n = —4. 
This confirms that in the relatively thin slice C < A both n = —11/3 and n = —4 cases, irrespectively of n have slope 
and amplitude determined by velocity fluctuations, as comparison with bottom left panel reveals. The slope is equal to 
—8/3, if the slice is thin i.e. (£/A)|(K|A) m / 2 < 1, and tends to —10/3 for thicker slices. Only slices much thicker than 
the velocity scale, i.e. C >• A, reveal the underlying density spectrum. 



fact that changes of underlying density spectrum do not alter the observed spectrum for |Kr c | ^> 1 and 
C/r c < 1 confirms our claim (see Figure 2). 

Now it is time to address the question whether our results depend on the cut-off size r c . For this 
purpose we plot the results of calculations for very different values of r c j A in Figure 3. It is easy to see 
that only marginal dependence of the spectral slope on r c exists for sufficiently large values of |K|. 

All in all, in the case of long-wave dominated density statistics both thin and thick slices provide 
velocity spectral index. To determine the density spectral index one should use very thick slices with 
C > A. Only then the emissivity spectrum begins to reflect the index of the actual HI spectrum 
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of density (i.e. the density spectrum in Galactic coordinates). The regimes of long and short-wave 
dominated spectra may be distinguished by varying the thickness of velocity slices. 




1 10 10 2 10 3 10 4 



KA 

Fig. 3. — Dependence of the two dimensional power spectrum on the structure function cutoff r c . The results for four 
values of the slice thickness L = C/X and three values of r c are shown. In shortwave regime |K|r c > 1 2D spectrum is 
completely insensitive to r c for thin slice and exhibit only minor dependence on r c if slice is thick. The latter dependence 
reflects the fact that thick asymptotics is achieved a bit later for smaller value of r c . 



4.3. Warm and Cold HI 

When we deal with cold (T m lOOfT) HI, the turbulent velocity is supersonic, while for warm (T m 
6000 K) HI the turbulent velocity is comparable with the thermal velocity. Our studies above dealt with 
supersonic turbulence and the effect of the thermal velocities was shown to be equivalent to increasing 
the thickness of the slice from ~ 5V to ~ 2 * (5v 2 /Q + 2-yf,) 1 / 2 (see section 3, we use W" = 1/6 
assuming step-like experimental window). As a consequence there is no point of getting experimental 
velocity slices narrower than Vl2vT,coid ~ 2.6km/ s. For warm HI the turbulence corresponds to the 
Mach number of the order of unity and the minimal thickness of the slice is ^ 17 km/s. What are the 
observational consequences of this? 

First of all, we may attempt to estimate the actual temperature of HI. We calculated how spectra 
should change with the thickness of the velocity slice. When 5v ^> vt the change of the 5v should 
entail the corresponding change in the effective slice thickness. However, as 5v approaches vt a further 
increase of the instrument velocity resolution (i.e. decreasing 5v ) will not result in the change of the 
power slope. In the medium with a few HI phases having different temperatures Ti the transitions will 
happen for a number of 5jV = 3.5^- Modulations of power-law spectral index as a function of 8V 
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should reflect the distribution of temperatures of HI components along the line of sight. We shall discuss 
the particular technique of probing HI temperature distribution elsewhere. 

For the generally accepted picture of the mixture of warm and cold components the power spectrum 
can be presented as the sum of a three: 

P 2 (K) = P 2 

(cold) 

(K) + P 2 

(warm) 

(K) + P 2 

(warm— cold) 

(K) (51) 

where the part P 2 ( W arm-coid)(K) arises from correlations between warm and cold components. Note that 
the minimal effective thickness of the slice for the warm component is ~ 20 km/s. When Sv approaches 
this thickness we expect to see variations of the spectral index if the warm gas constitutes a substantial 
fraction of the total HI. As mentioned above, this potentially may be used to determine the relative 
abundance of the components. 

Assume, for the sake of simplicity, that the mass fraction of the cold and warm media are comparable 
and the same is true for the amplitudes of the three dimensional spectra P s ( war m) and P s ( co id)- When 
the experimental slice width SV is smaller than Vl2vT( W arm) the the ratio of cold to warm contributions 
Pa (warm) / P a (cold) as a function of wavenumber K will be constant of order unity at long wavelength 
K < {C 2 /3.5v T{warm) )^ m (see criterion fl36|)), then fall off as K~ 2 ' m until K « {C 2 /5Vfl m and be 
constant again but suppressed as ~ SV/3.5vT( wa rm) at shorter wavelength (if 5V < 3.5t>T( co zd)), until the 
maximum suppression ~ VT( co id) / VT(warm) is reached) Q. These three regimes correspond to wavelengths 
for which effective slice is, first, thin for both cold and warm component, then still thin for cold, but 
already thick for the warm one and, finally, thick for both contributions. Addressing the actual Galactic 
HI data in Fig. 4 we may see that the slicing is frequently thin for cold HI and thick for warm HI. As 
SV increases above thermal velocity of warm component, shortwave relative suppression of contribution 
from warm HI disappears. 

In reality P s ( war m) and P s (coid) m ay have different amplitudes and even different spectral slopes. More- 
over, the mass fraction in those component may vary. The variations P 2 with the slice thickness may 
allow to determine the characteristics of the turbulence in the multiphase media. 

The P 2 (warm-coid) contribution arises from the correlations between fluctuations in warm and cold 
media. The density fluctuations are spatially separated in two media and therefore their correlation is 
likely to be negligible. We may not be so sure about the velocity fluctuations as cold and warm gas may 
participate in a coherent motions. It may be shown, however, that the difference between VT( W arm) and 
VT(coid) will suppress a possible contribution from P 2 (warm-cold) in a manner similar to the suppression 

of Pi (warm) • 

5. OBSERVED SPECTRA 

Revisiting L95 

In L95 the relations between 2D statistics available via observations and the underlying 3D statistics 
were obtained on the assumption that the observed statistics is determined by density fluctuations. 
Below we use our results on velocity mapping to interpret power spectra of intensity obtained by Green 
(1993) in terms of HI density and velocity fluctuations. 

Green's observations of the HI emission were accomplished with the Synthesis Telescope of the Do- 
minion Radio Astrophysical Observatory (DRAO) towards I = 140°, 6 = 0° (03 h 03 m 23 s , +58°06'20', 
epoch 1950.0) and they revealed a power law spectrum of 2D intensity. This spectrum is proportional 
to P 2 (K)\c and its interpretation depends on whether the slicing is thick or thin. 

To answer this question we should estimate both A and |K|. To estimate A we use a crude estimate of 
/, namely, we describe ordered motion of atomic hydrogen in the outer parts of the Galaxy by a simple 
flat rotation curve without any distorting radial motions. In this case the parameter / which describes 
line-of-sight projection of the shear in a rotational flow is 

/ « _ J ^ (52) 

2A (a + sinZ) 2 ,/! - (a + sin/) 2 



22 0ur earlier statement of the possible change of the spectral index can be understood as through this suppression. 
Indeed, for SV/3.5vT( W arm) 1 the contribution of the warm media is negligible, while it increases with SV. 
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where A = 14 km/s/kpc is the Oort's constant and a = V z /Vq, where Vo ~ 220 km s _1 and V z is the 
relative velocity of HI parcel. This relative velocity varies from at z = to —Vq sin(Z) for the gas 
parcel at z = oo. Correspondingly, the parameter /, which is negative, grows in magnitude from the 
familiar l/v4sin(2Z) in the Sun vicinity to larger negative values for distant HI regions. In particular, 
Green's data covers the velocities V z up to ~ —100 km/s, thus / ranges from —72.5 pc/(km/s) to —460 
pc/ (km/s). 

If we assume that velocity variations at the scale 30 pc amount to 10 km/s and arise from the 
Kolmogorov turbulence, the structure functions of velocity are 

/ \ 2/3 

D LL (r)^100f^— J kmV 2 , (53) 

where the cut-off of the function at large scales is disregarded on the assumption that turbulence is 
being studied at the scales much less than r c . Thus A given by eq. (|TJ) is ~ 3.5 kpc for the closest slice 
and ~ 56 kpc for the most distant slice. 

The width of the interferometer channels combined to give a single data point in Green's dataset is 
5V = 5.94 km/s. The slice thickness in parsec is £ ~ 5Vf pc, and varies from « 600 pc for the closest 
slices to ~ 2200 pc for the distant ones^. The wavenumber of transition from thin to think slice given 
by eq. (|36]) is equal 0.16pc _1 . 

Figuref| shows the range of the scales involved (see Green (1993)). The smallest |K| span from 
~ 1/3 pc -1 for the closest slices to 1/200 pc -1 for the distant onesQ, and therefore |K|A 1. The 
range of values corresponding to |K| _1 is shown by the darkened region in the plot. If we deal with cold 
HI, most of them, except the nearest slices correspond to the regime where slice is thin. For warm HI 
the border line is different and most of the scales correspond to very thick slice. 

We have shown earlier that the contribution of the warm component is suppressed unless the slices 
are thin for both components. The latter regime corresponds to the parameter space below the line 
"cold HI" in Fig. 4. For such close slices the actual geometry of the diverging lines of sight must 
be accounted for (see Lazarian 1994b) and we consider the resulting complications elsewhere. If the 
intensity contributions of the warm and cold HI are comparable, it is possible to show (see section 4.3) 
that the cold HI contribution will dominate the spectra above the "cold HI" line. 

Our analysis of Figure f| and Figures 1 and 2 shows that if the spectrum of density is shallow (i.e. 
n > —3), the observations by Green (1993) reveal the spectrum with index n + m/2. For m = 2/3 
the spectrum of emissivity obtained by Green (1993), namely the emissivity with the index ~ —2.7, 
corresponds to n ~ —3. This conclusions is not sensitive to the values of density correlation scale ro 
(see Eq. fl3"0"D), which we estimated to be of order of 20 pc. 

If, however, the density spectrum is steep (i.e. n < —3), the fluctuations of 21 cm intensity observed 
by Green (1993) can arise from velocity fluctuations. In this case the spectral index is —3 + m/2. For 
m = 2/3 one gets the slope —8/3 ~ —2.7 which is exactly what is observed. 

In the case of the Galaxy, the steepening of the spectrum depends on the scale |K| _1 as shown in 
Figure 5. This figure should be not interpreted literally for large C, however. First of all, when the 
thickness of the region becomes comparable to the distance to it, one should account for the divergence 
of lines of sight (Lazarian 1994b). Moreover, our model that assumes infinitely extended media does 
not hold and effects of finite size of the emitting region become important when L becomes comparable 
with the size of the Galactic disk (see Appendix E). 

SMC 

Small Magelanic Cloud (SMC) is rather irregular galaxy and the Earth observer looks at it from outside. 

23 Note, that the cut-off due to thermal velocity (see section 2.1) in Warm Neutral Medium (see table of idealized phases 
in Draine & Lazarian (1999) ) is ~ 6 km/s. If the WNM constitutes the dominant fraction of the neutral phase (Dickey 
1995) then the velocity resolution above is optimal and no further decrease in 5V will result in getting new information. 
However, if close to Galactic plane Cold Neutral Media constitutes a substantial portion of mass, the increase of velocity 
resolution up to 1 km/s is desirable. 

24 We stress that our analysis should be modified for the largest scales which are larger that r c . 
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-V (km/s) 

Fig. 4. — The variations of geometric scales with the sampling velocity are shown. The upper curve corresponds to 
the variations of the correlation scale A in the velocity space. The middle curve corresponds to the variations of the slice 
thickness C. The darkened area in the Figure depicts to the range of the turbulence scales under study in Green (1993). 
The solid curve within the darkened area corresponds to the interferometric measurements with the baseline 21m. The 
lower horizontal line denotes value of |K| _1 which separates thin (above) and thick slice regimes for cold HI and the 
upper line denotes a similar value for the warm HI. 
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Fig. 5. — The variations of the spectral index n of the 21 cm intensity spectrum with the thickness of the slice C for 
the turbulence with Kolmogorov velocity spectrum and long-wave dominated density spectrum. To provide the feelings 
of the scales involved we expressed our results in parsecs. 



For its study the modification of our technique described in the Appendix E is applicable. As the result, 
the size of the SMC plays the role of the parameter A and with this modification our earlier formulae 
are applicable. 

An advantage of using SMC is that the effects related to the divergence of lines of sight are negligible 
for the study and this simplifies the analysis for thick slices. Data on SMC in Stanimirovic et al. 
(1999), Stanimirovic et al. (2000) and Stanimirovic (private communication) shows that for thin slices 
the slope is ~ —2.8 but it steepens up to ~ —3.4 as the integration over the whole emission line is 
performed. Note, that in the latter case the intensity fluctuations arise from density fluctuations only 
and the spectral index of P 2 corresponds to the 3D spectral index of underlying density turbulence. This 
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means that the density spectrum has an index —3.4, which is close to the Kolmogorov —11/3 index. As 
this index is less than —3, the fluctuations should be "velocity-dominant" and the corresponding index 
should be — 3 + m/2, provided that the slice is "thin". The measured spectrum is a bit steeper than 
in the Galaxy. Is this difference meaningful? Does it mean that the spectrum of turbulence deviates 
from the Kolmogorov law? Is this a reflection of a transition to thicker slicing? A detailed study of the 
SMC data should answer these sort of questions. At the current stage it looks to us that there exist a 
rough agreement between the data and the assumption that both density and velocity spectra have the 
Kolmogorov index of —11/3. 

An additional interest to SMC data stems from the fact that variations of the power law may shed light 
on the nature of turbulence in warm HI. We discussed above that the fluctuations in sufficiently thin 
slices are dominated by the cold component. As the thickness of slice increases the relative contribution 
of the warm component increases. The fact that the velocity-integrated data provides the index close to 
the Kolmogorov one either indicate that the fluctuations of warm HI density are insubstantial or that 
they have the Kolmogorov index. 

Dust and Molecular Data 
Dust and HI are well correlated (see Boularger et al. 1996). Both Galactic and SMC Far Infrared emis- 
sion exhibits shallow spectrum with the index —3 (Wall & Waller 1998, Waller et al 1998, Stanimirovic 
et al. 2000). However Stanimirovic et al. 2000 combined data at 100 /zm and 60 /im to show that the 
actual dust spectrum of dust column density in steeper. In fact the spectrum of dust column density 

was found has index « —3.6 which is in rough correspondence with the ~ —3.4 3.5 spectral index 

for HI (see above). Thus both HI and dust supports the idea that the spectrum of density is steep for 
diffuse media. 

Are there any evidence that that the interstellar density may be shallow? In the theory domain 
Henriksen (1994) derived a quite shallow spectrum for compressible turbulence, but we do not know 
whether this result is applicable to HI. Although molecular cloud statistics may be different from that 
of diffuse HI, comparison of the two is worth doing. 

Recently molecular cloud data was discussed in Stutzki et al. (1998). There for both 12 CO (data from 
Heithausen & Thaddeus (1990) and Falgarone et al (1998)) and 13 CO (data from Heyer et al. (1997)) 
transitions the spectrum of intensity was observed to have a power law index ~ —2.8. As the data is 
averaged over velocity, the fluctuations of intensity are due to density fluctuations and the spectrum 
of density should have the same slope as the intensity (i.e. n = 2.8, (see Eq. (p7|)), provided that the 
transitions are optically thin. The trouble is that the lines are optically thick and the interpretation 
above is not applicable. Further research should determine what the particular emissivity intex means 
in terms of underlying statistics of velocity^ and density. 

Strangely enough, the shallow spectrum of molecular density fits well in the picture of molecular cloud 
structure. Stutzki et al 1998 derive the relation [^between the power spectral index n, the spectral index 
e of clump mass-size relation M cx r e , and the index K of the clump mass spectrum dN/dM oc M H : 

n = (3 - N) e 

e = 2 corresponds to the Larson's relations (Larson 1992) and K in the range 1.6 to 1.8 corresponds to 
the analysis of clumps (Kramer et al. 1998). The trouble is that this relation was derived to be true 
for optically thin clouds, while the clouds in the analysis by Stutzki et al 1998 are optically thin makes 
one worried. Further research should concentrate on providing the adequate statistical description for 
optically thick tracers and comparing the statistics for thin and thick tracers. 

6. DISCUSSION 

25 Although the data is integrated over the emission line, the random velocity is important as self-adsorption depends 
on the gas velocity dispersion. 

26 This relation can be applied to individual structures in velocity slices provided that the index of the two dimensional 
spectrum of intensity fluctuations is used instead of n. In the velocity dominated case the masses would reflect not the 
masses of the real clumps but caustics produced via projection. 
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6.1. Our approach 

Velocity mapping In the paper thus far we address the problem of the velocity modification of HI 
statistics. Although the paper deals with HI only, the problem it addresses is quite general and therefore 
our results are applicable to any optically thin lines. We plan to address the effects of finite optical 
depth elsewhere. 

We showed that intensity fluctuations can arise both from velocity and density fluctuations and derived 
the 3D spectrum of intensity (eq. (0)) within velocity data cubes. A quick look at eq. (|T3| ) shows that 
the density and velocity enter the expression in a different manner and this eventually enables us to 
separate the velocity and density contributions. 

Statistics used: 2D Spectra 
Channel maps and the related 2D spectra constitute the central topic of this study. This is because the 
channel maps are readily available via single dish observations, while channel map spectra -P2(K) can 
be directly obtained with an interferometer. 

As we discussed in the previous section it is important to determine whether 21 cm intensity fluc- 
tuations are due to density or velocity fluctuations. If the density spectrum is long-wave dominated, 
our results in Figure 2 indicate that HI intensity spectrum does not reflect the underlying HI density 
statistics. Therefore the spectrum of the density fluctuations may have any value n < —3, but this will 
not alter the intensity fluctuations unless a substantial portion of the 21 cm line is integrated over and 
the slice become "thick" in accordance with the criterion (j36l). 

Additional statistics: ID and 3D Spectra 
In view of the problems we face with the interpretation of present HI results it is important to make 
use of other statistical tools. Indeed, 2D statistics of channel maps is not the only statistics that is 
available. Using data cubes one can obtain the 3D spectrum P s which properties are described in the 
Appendix B. Alternatively, an additional insight can be obtained if ID spectra are formed using velocity 
cubes. We discuss these new statistical tools in Appendix F and show that ID spectra is complementary 
to the tools we employed in the main text. In fact, using this statistics one does not need to change 
the velocity slice thickness to distinguish between long and short-wave dominated density spectra. If 
the density field is long- wave dominated we expect to see the power index — 2/m, while if the density 
spectrum is short-wave dominated we expect to get the index 2(n + 2)/m. This test can be performed 
with the existing 21-cm data cubes. For n = —3 the predictions for the two cases coincide. Naturally 
the spectra are not independent, but they are produced using different procedures of data handling and 
therefore provide an additional test. 

6.2. Problems of Interpretation 

"Big Power Law" 

The application of our technique to the available data (see section 5) testifies that the density spectrum 
of HI in SMC has a power-law index close to the Kolmogorov value of —11/3 over the scales from 40 pc 
to 4 kpc. Assuming that density statistics in our Galaxy is steep (as in SMC) we also get that the 
velocity has an index of ~ —11/3 from 10 to 100 pc. The proximity of indexes for both velocity and 
density spectra to the Kolmogorov values is both remarkable and unexpected. First of all, there are good 
reasons why we should not a priori expect to find Kolmogorov-type turbulence in interstellar medium 
(Lazarian 1999a). Moreover, Kolmogorov interpretation would entail problems with understanding of 
the nature of energy injection. We also noted earlier that SMC possibly shows a density spectrum with 
only marginal velocity contribution, which would mean that the structures we observed are static. The 
last point requires a further investigation, however. 

It is very strange, but Astrophysical turbulence whenever it is measured seem to exhibit Kolmogorov 
power slope. The existence of the "Big Power Law" (see Armstrong et al 1995, Spangler 1999) is 
one of the great Astrophysical mysteries. Mesurements of turbulence in ionized media established a 
spectrum with —11/3 index at the scales from a few AU to a fraction of a parsec. Measurements at 
larger scales are much less reliable (Cordes 1999) and although compatible with the —11/3 index allow 
other interpretations. Contrary to our intentions (see L95, Lazarian 1999a) we contributed to making 
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the problem even less tolerable. It looks now that the "Big Power Law" spreads up to scales of several 
kpc. 

MHD turbulence 

Ignoring for the moment the issue of energy injection scale and associated problems we may ask whether 
we expect to see the turbulence with the index close to the Kolmogorov one when the fluid is strongly 
magnetized. 

The Kolmogorov-type spectrum of plasma density fluctuations observed via radio scintillations and 
scattering (see Armstrong et al 1995 and references therein) has been interpreted recently as the conse- 
quence of a new type of MHD cascade by Goldreich & Sridhar (1995). The Goldreich-Shridhar model 
of turbulence^ differs considerably from the Kraichnan one (Iroshnikov 1963, Kraichnan 1964). It 
accounts for the fact that hydrodynamic motions can easily mix up magnetic field lines in the plane 
perpendicular to the direction of the mean field. Such motions provide eddies elongated along the field 
direction and the velocity spectrum close to the Kolmogorov one. 

The Goldreich-Shridhar turbulence is anisotropic with eddies stretched along magnetic field direction. 
The wavevector component parallel to magnetic field k\\ scales as k 2 ^ 3 , where k± is a wavevector compo- 
nent perpendicular to the magnetic field. Thus the degree of anisotropy increases with the decrease of 
the scale. The fact that during the observations various scales contribute to the correlation functions of 
intensity (see discussion in L95) and the contribution of the largest scales dominates the signal mitigate 
the anisotropy of the intensity statistics. The component of magnetic field parallel to the line of sight 
decreases anisotropy as well. Indeed, it is easy to see that only the component of magnetic field that 
is in the plane of the sky will create the anisotropy, while the component of magnetic field parallel to 
the line of sight "mix up" parallel and perpendicular direction^]. Therefore the anisotropy of the 21-cm 
statistics may not be easy observable (see Green 1994). In our future papers we plan to discuss the 
observational tests of the Goldreich-Shridhar picture for HI. A simultaneous use of other techniques of 
turbulence study (see Lazarian 1992, Lazarian 1993) should allow to have a more comprehensive picture 
of MHD turbulence. 

Incompressibility is assumed within the Goldreich-Shridhar picture of turbulence. In actual ISM 
turbulence the compressibility may be important in view of the fact that the media may be pliable 
to compressive forces (soft equation of state). It is unclear how the picture of MHD cascade will be 
changed in this case. 

When magnetic field is intermittent, one may expect the thermal instability to be more efficient in 
the regions with a smaller value of magnetic field. Indeed, if varios regions of magnetized gas are at 
pressure equilibrium, the gas pressure within the regions with lower magnetic field will be higher. This 
shall increase the density of gas and may initiate thermal instability. Indeed, the heating of gas is 
proportional to the gas density, while cooling is proportional to the (density) 2 . As the result dense gas 
will tend to cool down further increasing its density (Spitzer, 1978) making regions with lower magnetic 
field colder and denser. This would qualitatively correspond to the picture of magnetic measurements 
in HI by Carl Heiles (private communication). We will quantify the picture and its implications on the 
small scale structure of HI elsewhere. 

Superbubbles and Supershells 
One may speculate that the observed statistics is the result of the superbubbles and supershells. Indeed, 
potentially superbubbles may inject energy on the kpc scale. The issue in this case is why we do not 
see "knees" in the spectra at the intermediate scales when singular shells inject energy. Does this mean 
that the supershells should carry an overwhelming amount of energy or the available resolution is not 
sufficient to see peculiarities related to more localized energy injection? 

In general, energy injection scale is critical for understanding interstellar turbulence. For decades 
it was generally believed that MHD turbulence damps slowly (Iroshnikov 1963, Kraichnan 1964). In 
this situation supersonically broadened linewidths of molecular clouds could be explained as locally 

27 A qualitative discussion of the model and the role of reconnection for the cascade can be found in Lazarian & Vishniac 
(1999). 

28 No anisotropy is expected for the intensity correlations if the total magnetic field is parallel to the line of sight. 
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generated, e.g. due to a cloud collapse. However, new understanding of MHD turbulence (Goldreich 
& Shridhar 1995) and numerical simulations (e.g. Vazquez-Semadeni 1999) suggest that the turbulent 
damping in magnetized media takes place on the scale of an eddyQ turnover (similar to the case of 
unmagnetized media). Therefore to allow starless cores to have linewidths similar to those of star- 
forming cores (Benson & Myers 1989) the energy injection should happen on the scale much larger than 
the cloud size (for larger eddies the timescale of turnover is larger). Nevertheless, scales of the order of 
several kpc look excessively large. 

Energy injection to molecular clouds may arise directly from shocks (Scalo Sz Kornreich 1999), and 
this also can be classified as injection from large scales. The spectrum of shocks is steep (~ k~ A ) and 
does not correspond to the observations. 

One have to accept that power spectra may be rather rough statistical tools and the index —11/3 
may have nothing to do with a cascade. For instance, such a density spectrum may emerge in a media 
with HI shells having a distribution P(x) ~ x s ^ 3 dx. However, the particular power law does not seem 
to have particular physical explanation and moreover, small shells are expected to have larger velocities 
which does not correspond to the observed velocity spectrum. 

6.3. Model assumptions and future work 

This work is the first, as far as we know, that treats quantitatively fluctuations arising from turbulent 
velocities. It is not surprising therefore that there are many issues to be studied in depth in the future. 

Density- Velocity Correlations 
Our derivations above assumed that the turbulent velocity and density are not correlated in Galactic 
coordinates. This seems to be true at large scales, but may fail on small scales where self gravity may 
contract HI structures^. Such a correlation may enhance the amplitude of fluctuations observed if 
the local divergence of the velocity field is negative in the regions of density enhancement (i.e there is 
systematic inflow of matter into overdense regions) and, vice versa, would decrease the amplitude of 
fluctuations if overdense regions are associated with outflows. However, interstellar medium is far from 
being quiescent even on the small scales and the sense of velocity-density correlation is unclear for the 
ISM turbulence. It is easy to see from the equation of continuity that the logarithmic time derivative 
of density equal to the divergence of the velocity field. This means, that the instantaneous density 
itself depends on the integral of the velocity divergence and in the complex medium with sequential 
expansions and collapses the importance of this correlation is difficult to estimate. In any case, we do 
not expect a priori any strong velocity-density correlations. 

To check whether the assumption about density-velocity correlation is important for our final results 
one needs either to do numerical simulations or to assume the statistics for joint density-velocity dis- 
tribution. In Appendix D we check our results assuming Lognormal distribution for random density^ 
and the velocity and logarithm of density are correlated Gaussian fields. 

To estimate the upper limit of the effect of density- velocity correlation we assume the maximal possible 
value of this correlation. Naturally, in all realistic situations we expect the value of this correlation to be 
much smaller than this maximal value. Our results for this worst possible regime can be summarized as 
follows. The criterion for the slice to be thick is increased at most by a factor of a few and all our findings 
for the thickF] slices (see Table 1) hold. The correspondence between our earlier results and "extremely 
correlated" results holds for the thin slices in the long-wave regime, which is favored by our analysis of 
data. For thin slices and short-wave dominated density field (shallow density spectrum) the correlations 
are shown to enhance the velocity effects. In fact, in the extreme regime of correlations it is velocity, 
rather than density that dominates the expected channel emissivity. However, we stress that we used 

29 It is shown in Lazarian & Vishniac (1999) that MHD turbulence allows equivalent descriptions via wave-wave, wave- 
eddy and eddy-eddy interactions. 

30 Whether this important in practical terms of HI studies is unclear, as a transition to higher densities entails a transition 
to molecular hydrogen. 

31 Lognormal distributions naturally arise in numerical simulations of isothermal gas (see Vazquez-Semadeni & Passot 
1999). 

32 This includes both long-wave and short-wave density field. 
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an upper limit for these correlations and for interstellar data sets we expect to observe density effects. 
Changes of the spectrum within the thin slice regime when the density is short-wave dominated can be 
a signature of strong velocity- density correlations. Using the thick slices it is possible to check whether 
the density is short-wave or long-wave dominated. A more detailed study of the regime of velocity- 
density correlations can be done using our expressions in Appendix E and substituting a different law of 
velocity-density correlations (rather than taking its maximal value), but this study is beyond the scope 
of the present paper. 

Anisotropy and Linear Expansion 
It is both challenging and important to determine the degree of the HI spectrum anisotropy and its 
variation from slice to slice. This information can provide an insight to the nature of HI turbulence and 
may be used as a diagnostic for the interstellar magnetic field. For instance, measuring the structure 
functions of HI intensity as a function of a positional angle in an analogy with the statistical treatment 
of synchrotron fluctuations (see Lazarian 1992) may reveal magnetic field direction in various portions 
of the sky, provided that the turbulence is indeed Goldreich-Shridhar type. So far, the attempts to 
measure anisotropy in HI are limited to the Green (1994) study, where no anisotropy was detected. 
Apparently a better analysis is needed. For the slices with high degree of anisotropy our statistical 
technique should be modified as suggested in L95. 

Another issue is related to the linear expansion (fj). If the slice is sufficiently thick the map becomes 
non-linear. Potentially, knowing a detailed Galactic rotation curve one can describe a non-linear map- 
ping. We found, however, that the linear mapping is OK for both thin and thick regime if we study 
sufficiently small scales (see discussion after eq. (1)). We also found that if L is smaller than the distance 
to the slice, it is possible to ignore the complication related to the radial nature of lines of sight. Thus 
the theory at its present stage is widely applicable. 

Applicability of the theory 
The theory presented in the paper is quite general. The mathematical machinery developed here can be 
used for studying turbulence in various emission/adsorption lines. For instance, velocity effects and line 
of sight integration are the problems that are also encountered in studies of turbulence inside molecular 
cloud. Unfortunately, there are additional complications related to studying molecular clouds. For 
instance, these studies face problems with samples being not statistically homogeneous (see Miesch 
& Bally 1994). Indeed, it is widely believed that for individual molecular complexes the statistics 
(especially at large separations) is dominated by regular gradients rather than the random component. 
To eliminate the inhomogeneous component Zurfleh (1967) or other types of spatial filtering (see Spicker 
& Feitzinger 1988a,b) can be used. However, Stutzki et al (1998) have shown that the statistics can be 
obtained without such filtering if the particular modification of wavelet analysis, namely, so-called Alan 
variances are used. For power-law statistics Stutzki et al. (1998) relate their statistical measures to the 
spectrum and correlation functions that are used in the present study. Therefore the translation of our 
results to the language of Alan- variables is trivial. 

What is less trivial is to account for the finite optical depth effects for particular emission lines. In 
treating this problem the adsorption by dust presents the least of the evils. We expect more complex 
dependences treating self adsorption within an emission line. The optical thickness in this case will 
depend on the velocity adsorption. 

Our present results are obtained for 21 cm emission transitions which intensity is proportional to 
the integrated density of HI along the line of sight. Other emissivities, e.g. those of H a lines are 
proportional to the squared density of species. It is possible to generalize our results for those transitions 
and contribute to the studies of ionized emitting media, e.g. of HII regions (O'Dell 1986, O'Dell & 
Castaneda 1987). 

Other generalizations are also straightforward. Our present research was motivated primarily by the 
necessity to interpret the existing interferometric data and therefore in the paper we mostly dealt with 
spectra. As the result, correlation functions (see Appendix C) serve auxiliary purposes in our study. 
With other sets of data, correlation functions may be more advantageous and our analysis of their 
properties in Appendix C is useful. Moreover, modifying our technique it is possible to study various 
moments of the line-shape (Lazarian 1992), e.g. velocity centroids (see Dickman 1985). In our future 
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paper we shall discuss all these possibilities. 
HI data and simulations 

In the paper we predicted how the power slope will change with the thickness of the velocity slice as 
a function of velocity and density statistics and the ratio of the cold and warm HI filling factor. It is 
important to test these our predictions against observations. The SMC data is an example of the data 
set that can be used for the purpose. 

It is also important to verify our ideas using simulations. Our results on velocity and density correlation 
were tested so far only for Lognormal distribution where we managed to get an analytical result. Using 
the statistics that arises from simulations is on the agenda. We would stress that the only 3D simulations 
are useful. Due to a limited dynamical range these simulations do not show power-law behavior and 
additional care is required for comparing these results with our predictions. 

Clouds and Filaments 

HI data cubes, exhibit a lot of small scale emissivity structure^. The question is what part of them 
is real, i.e. is associated with density enhancements in galactic coordinates and what part of them is 
produced by velocity fluctuations. A related question is whether the structures we see are produced 
dynamically, through forces, e.g. self-gravity, acting on the media or they may be produced statistically 
exhibiting the properties of random field. The second question was partially answered in Lazarian & 
Pogosyan (1997), where it was shown that density fluctuations with Gaussian distribution and power 
spectra result in filamentary structures. The structures become anisotropic and directed towards the 
observer when the velocity effects are accounted for. 

The issue of density enhancements produced by velocity fluctuations is closely related to the statistics 
of "clouds" observed in velocity space. The results on velocity mapping that we discussed earlier suggest 
that spectra of fluctuations observed in velocity slices are more shallow than the underlying spectra. 
This means more power on small scales or, in other words, more small scale structure ( "clouds" ) appears 
in the velocity slices due to velocity fluctuations. Understanding of the velocity mapping may help to 
solve the problem of discrepancy between counts of clouds observed in emission and extinction (compare 
Scalo & Lazarian 1996). 

The integrated velocity maps are sensitive to density fluctuations and may provide means of solving 
the problems of the statistics of actual density inhomogeneities. 

Beyond 2-point Statistics 
Is the power spectrum the best one can get? Power spectrum does not reflect the clustering properties 
of the medium. To find those characteristics one should use multi point statistics, of which the simplest 
would be the so-called bi-spectrum (Scoccimarro 1997). This type of statistics is sensitive to clustering 
properties of the media and we may hope to relate the bi-spectrum with the inflows observed in molecular 
clouds (see Myers & Lazarian 1998) and the processes of star formation. Other potentially useful tools, 
e.g. genus analysis are discussed in Lazarian (1999b). A choice of a particular tool depends on the sort 
of questions one attempts to answer. This paper is limited to making use of channel map spectra. 

All in all, we believe that the machinery developed in the present paper opens avenues for both 
theoretical and observational research and the appropriate studies should contribute to unraveling the 
mystery of ISM turbulence. At the same time a search for new statistical tools appropriate for diagnostic 
of ISM turbulence should continue. 

7. CONCLUSIONS 
The results of the present work can be summarized as follows: 

1. The emissivity spectrum arises both from velocity and density fluctuations. The emissivity statistics 
depends on the thickness of the velocity slice. The minimal value of the effective thickness depends on 
the temperature of the gas. We found limiting regimes that we termed "thick" and "thin" velocity slices 
(for criteria see Eq. fl35D ) and obtained asymptotics for these regimes. Thermal velocity acts to increase 
the effective thickness of the slice without increasing the amplitude of the signal. For a mixture of cold 
and warm HI the slice may be thin for cold phase and thick for the warm phase. In this case the the 

33 It was noticed by Langer, Wilson & Anderson (1993) that more structure is seen in spectral line data cubes than in 
the integrated intensity maps. 
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expected relative contribution of warm medium is suppressed. 

2. Two dimensional spectrum of 21 cm intensity (that can be measured by interferometer) corresponds 
to the integral of the underlying 3D spectrum if the slicing is thin while in the thick slice regime the 
result depends on whether the density is short- or long-wave dominated (see below for the criterium 
definition). The 3D spectrum is anisotropic and can be expressed through the velocity and density 
statistics. 

3. For Galactic data we find that the slice is thick if the two dimensional wavenumber |K| > 
(6V 2 /C)^ m « 0.16 per 1 , which means that the present day radiointerferometers provide a thin slicing 
of Galactic disk data. Density and velocity enter our formulae in a non-symmetric fashion. Therefore 
the slope of the HI emissivity spectrum changes with the thickness of the slice and this allows to find 
the slopes of the spectra of velocity and density fluctuations separately. 

4. For thin slices the velocity fluctuations make spectra of emissivity more shallow creating a lot of 
structures in PPV space that can be erroneously identified as clouds. For long-wave dominated density 
most of the small scale structures are due to velocity fluctuations both in thin and thick slice regimes. 

5. If the overdensity is dominated by small scale inhomogeneities (i.e. the correlation functions scale 
as 1/r 7 , 7 > 0), the statistics of the channel maps is dominated by density fluctuations. The HI 
observational data seem not to support shallow density spectrum, however. 

6. If the overdensity is dominated by large scale fluctuations (i.e. the structure functions scale as 
1/r 7 , 7 < 0), the statistics of the channel maps is dominated by velocity fluctuations. In this case, to 
measure density one should use very thick slicing £ ^> A, where A is given by Eq. (ITBp. In thin slices 
one will find universal 2D spectral slope ~ —8/3 if velocity field is described by the Kolmogorov index 
m = 2/3. Data on the galactic HI statistics and that of SMC roughly agree with this spectrum. 

7. The two cases above (i.e. 7 > and 7 < 0) can be distinguished by varying the thickness of 
the velocity slice. The analysis of the observational data suggests that the HI density spectrum has 
an index close to the Kolmogorov one, namely, —11/3 and the statistics of the velocity field is also 
Kolmogorov-type. Thus one may argue that the observed emissivity spectrum might be a signature of 
the Goldreich-Shridhar (1995) MHD turbulence. The developed mathematical machinery is applicable 
to studies of turbulence both in Galactic disc and individual clouds (see Appendix E) . Various emission 
lines can be used and the technique can be modified to account for finite optical depth effects. 
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APPENDIX 



Relevant to Table [B^ 

3D power spectrum 

a n = 
b m = 



LIST OF NUMERICAL CONSTANTS 



27r 2 r(-n)/ 



r( 3±n 



r[i=* 



§ ttT(5 + m)/r (1 - m/2) T (3 + m/2) 



m/n 


-1 -5/3 -2 


-3 


1/2 
2/3 
1 


74.8279 434.657 1361.75 
35.4366 88.5916 172.584 
19.1096 23.0004 29.4949 


84244.6 
2205.93 
92.5001 



eq. dETpj) 



eq. ([B7|) 
eq. QBig) 
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where the numbers at the right hand sight correspond to the formulae in Appendix B where the constants 
appear first. 



Relevant to Table [I] and Table 

2D power spectrum, 

B rn = (13/57r) b rn 

s nm = 2 - n -( m+1 )/ 2 7rr (-3s±2i) / |r ( 3+n ^ m/2 ) 



Correlation function, 



2 ( 7+ m-l)/m r ( 2 ^~ 2 ) ^1 



c im = r (attest) / 2 lr (^± m 



ASYMPTOTICS AND SPECIAL REGIMES FOR 3D POWER SPECTRUM 
We discuss asymptotics and introduce useful approximations of the integrals involved in computation 



of the 3D power spectrum P s (k) = P P (k) + P v (k). If the density correlation function is given by eq. (|30D 
Introducing new dimensionless variables 

y = kr (Bl) 
x = cos 9 

<P = (k\) 2 - m n 2 = X 2 - m k 2 z k~ m , 
we can present the 3D spectrum P s (k) (see eq. (|T6|)) as 



P s {k)=P v (k) + P p (k) , (B2) 
P v {k) = k^l^,4>) , 
P p (k) = F- 3 rJJ_ n (/i,0) . 

lab, 4>) = 2n y^dy £ dxJ (yj(l - - a*)) cos^e"^^ 1 -^/ 4 , (B3) 

where a = 3 or — n and Jo is the Bessel function of the zeroth order. We have to calculate P v and P p 
in another combination when the density correlation function is given by eq. ( |4?D 
It is most illuminating to start by expanding the plane wave in the eq. ([13]) into multipoles 

47r^i I ii(A;r)y lm (k)^(f) , (B4) 



e ikr 



Im 



where j\(kr) is a spherical Bessel function and Y\ m are spherical harmonics. After performing trivial 
integration over the polar angle in eq. ([L3j) and making the dy integration the innermost one we get 

l a {k) = 8vr 2 J2(i})Y l0 (ii) £ dxY*(x) (^°° y^dyme^ 2 ^ 1 -^'^ (B5) 

In this Appendix we consider two asymptotics, one "long-wave" corresponding to small <fi and the 
"short-wave" corresponding to large (p. In the long- wave limit <f> (see eq. flB2|) ) tends to zero and the 
integration is trivial. In this case 1 a is just a constant which coincides with the value of X a in the 
absence of random velocity 

r A°)° h ' r( J(rte- ' < B6 > 



2 J V 2 
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The result is the same if one considers the wave modes orthogonal to the line of sight (/x = 0) which are 
unaffected by velocity mapping Z a (0, 0) = X Q (/x, 0). 



For the term P p arising from density perturbations 



a. 



-n 



< 3 and 2L n (/i, 0) is finite. For the 



velocity term P v , a = 3 and eq. ( |F10|) gives zero value. The 5-function at k = which is contained in P v 
and describes the mean density in velocity space has been implicitly omitted by the choice of a variable 
y. To get more informative result, we expand eq. (|B5|) up to terms with to the first power, and after 
integration over angles obtain 



Z 3 (>,0 -> 0) « b, n 



3 /l 



where 



m 9 

71"—^- 

12 r i 



1 + 5U-^ 



r(5 + m) 



m/2)r(3 + m/2) 



(B7) 



(B8) 



Note, when we expand the exponent in powers of each next term of the expansion subsequently gives 
rise to higher multipoles in I series. Namely, the zeroth-order term in leads only to a monopole 
contribution, the first-order terms add quadrupole, the second-order add an octupole etc. 

Let us derive the short-wave asymptotics for eq. ([Bo]). We note that the function ji(y) is peaked, 
roughly, at y — 1+ 1/2, so if the exponential cut-off occurs at lower y, the corresponding Z-term can be 
neglected. Thus only low multipoles remain in the sum for large 0. Indeed, for a given high multipoles 

-, \— 1/m 

I + 2 > ( gi ye vanishing contribution into eq.([B5|). Here we use an estimate for the average 

value of the angle < x 2 >= 1/3 to get (2 + m - m < x 2 >)/4 = ^±p. 

Short-wave (i.e. high 0) asymptotics follows when only monopole is left in the multipole expansion 
(|B~5|). This is valid, formally, when I = 2 quadrupole is cut out (since the integrals of a dipole is zero). 



The criterion for this is d> > 



3+m 



(!) 



1. Therefore 



T a (fi, > 1) « 2vr / y^dyUy) / dxe~^ 2+m ^ 2 ^ 
Jo J-i 

For very short waves, i.e. ^> 1 one can put jo(y) = 1 in eq.([B9|) to obtain the limit 



oo 



8tt / 4 ' 



2 d {vV^S 



D(y) is the Dawson integral. The factor in brackets, s nm , is tabulated in Appendix A. 



(B9) 



(BIO) 



We found that eq. ( |B9| ) well approximates expression ( |B5| ) for all 0. Indeed, both long-wave asymp- 
totics eq. ( |B6|) and the short-wave asymptotics eq. ( P39|) are obtained with only I = term in eq. (|B5| ). 
We checked numerically that eq.(p9|) provides an accurate approximation for X_ n , n ^ —3, for all values 
of /j, and 0. This is a very welcomed simplification, since eq. ( |B9| ) depends on rather than on k and 
H separately. For a = 3 at small the zeroth-order /^-independent term ( |B6| ) is equal to zero and the 
in linear term of X3 at low has the /i-dependent amplitude (see p7|) . In this case Eq. (|B9|) provides 
an angle-average estimate. 

The deviation of the expression ( |B9| ) from ( |B5| ) is less than 3% for a = 1. In fact, we found the 
following formula to be even more accurate numerically while still depending only on variable 



£*(!,■ 



2tt J y a ~ l dy J dxcos(yx)e~ 



■4>y™(2+ m (l-x 2 ))/4: 



(Bll) 



This expression is exact for the modes parallel to the line of sight (/i = 1) with = (k\) 2 ~ m . Being 
extended to general = (k\) 2 ~ m fi 2 it gives an excellent approximation to the full result. 

Above we discussed the asymptotics of the integral ( |B3[ ) that can be found analytically. The Table |B"2 
summarizes our results for the power spectrum in the velocity space. In velocity space the spectrum 
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(k z X) 2 < (kX) m 


(k z X) 2 > (kX) m 


A- 3 P P (K,^) : 


a n (r /X) n+3 


■ (kX) n 


Snm(ro/\) n+3 ■ {k z X) 2n ' m 


X- 3 P v (K,k t ): b m 


i + § (1/2 - ki/ej 


■ (k z X) 2 /(kX) 3+m 


s. 3m ■ (k z X)~ 6/m 



Table B2 

Asymptotical behavior of the 3D spectrum in velocity space. Numerical constants a n ,b m 

AND S nm ARE GIVEN IN APPENDIX A. 



is anisotropic, depending separately on 2D wavevector K, perpendicular to the line of site and k z 
component, parallel to the line of site, k 2 = |K| 2 + k 2 . 

The first ( "long- wave" ) asymptotics is achieved for waves longer than a velocity correlation scale A 
or for the waves transverse to the line of sight, which are unaffected by velocity mapping. The second 
( "short-wave" ) asymptotics is reached for waves that are shorter than the turbulence length and have 
non-negligible component in z-direction. 

In Fig. B6 and Fig. B7 we show results of numerical calculations for the total spectrum in velocity 
space P s = P p + P v obtained for r = A and r = 0.02A respectively. In the former case the transition for 
the short-wave asymptotics happens exactly at the same scale as density correlation scale, while in the 
latter, more realistic case the velocity scale significantly exceeds density correlation scale. Naturally, 
when density in Galactic coordinates is less correlated the caustics caused by velocity fluctuations 
become important. The comparison of Fig. 5 and Fig. 6 confirms this tendency. The results for 
P v and P p constitute the basis for our computation of 2D spectra in the main text and ID spectra in 
Appendix F. 

3D CORRELATION FUNCTIONS 

The two-point correlation function is an alternative to a power-spectrum description of the density 
field. In the velocity space the correlation function £ S (R, z) is obtained by the inverse Fourier trans- 
formation of the spectrum P s (k), £ S (R, z) = (2tt)~ 3 J d 3 ke -ik " rs P s (k). Here, as in the main text, we 
use capital letters to denote vectors perpendicualar to the line of sight and therefore the total vector 
in the velocity space is r s = (R, z s ) and R = |R|. With an anisotoropic P S (K, k z ) given by (|T3[), the 
integration over K in eq. (^l|) leads to delta function, which can be used to eliminate the integral over 
R. The integral over k z in eq. (|2lD can be done analytically to obtain 



UR^ S ) = Jdz£{R,z) (2nX 2 D z (R,z)) 1 exp 



2X 2 D z (R,z) 



(CI) 



Here D z (R,z) describes a velocity structure tensor, i.e. D z (R,z) = (r/X) m 1 + y(l — z 2 /r 2 ) in the 

case of solenoidal turbulence, and r 2 = R 2 + z 2 . In the absence of the random velocity field A —>■ the 
kernel in the brackets transforms to a delta function 5(z s — z), and, as one expects l; s (R,z s ) — > £(r). 
Similarly to the split of power spectrum, the correlation function can be split in two parts £ s = 
where unity describes mean density in velocity space, £ v comes from the random velocity map of the 
mean density term and corresponds to the first integral in (CI) minus unity, while £ p arises from the 
transformation of the underlying correlations of density perturbations £(?"). It is useful to introduce 
dimensionless variables f, R, z which measure distances in terms of the velocity correlation scale A, i.e 
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Fig. B6. — The quantity P(k)k 3 for various values of m and n = 7 — 3; ro = A. The short wave part of the spectrum 
is the function of the variable = (k z X) 2 /(kX) m only. Although the longwave spectrum depends separately on k z and 
k, we have shown the total spectrum in (f> variable, assuming the modes to be along the line of sight k z ~ k. In each 
plot the heavy line shows the total spectrum P s = P p + P v , while two light lines show P p and P v terms separately. The 
term originated from the underlying density perturbations P p dominates at large wavenumbers. 



r/A, etc. Then, with 7 = 3 + n denoting a slope of the correlation function £(r) 

z a ) = -1 + Jdz G(R, z s , z) (C2) 

£ P (R, z s ) = (r /\y J d~z (R 2 + -z 2 )^ G(R, ~z s , ~z) (C3) 

G(R, z s , z) = (2nD z (R, ~z))~ l exp f- ^'-fH (C4) 



Two specific cases are of particular interest, namely, 2D correlation function in a slice of fixed velocity 
C, S (R, 0) and ID correlation along the line of sight £ s (0, z s ). The first one corresponds to 2D correlation 



function measured in a thin slice. Asymptotical behavior of £ a (i£, z s ) is given in Table |C3 . 

The asymptotical formulae for £ p are valid for 7 + m/2 > 1. As it should be, the correlation function 
is unchanged by velocity map at large separations, having the same slope 7 as the underlying £(r). 
If 7 > 1, the slope of correlation function at small separations is changed, becoming more shallow 
—7 + (1 — m/2) for correlation in the slice of fixed velocity, but steeper —7 — ^ 2 ~^f > (7 — 1) for the 
correlation along the line of sight, if 7 > 1. Interestingly, for 7 = 1 there is no change of correlation 
slope along the line of sight. However the amplitude of correlation at small distances is enhanced by 
the factor 1/m. 
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(r/A) -> 


(r/A) — ► oo 


(A/r )^,(i2,0): 


C 7m • (i?/A)~ 7+(1 ~ m//2) 






C-/m ' \%8/ A) 


(z s /A) -7 




mj (2 — m) 


(m 2 /4)-(R/\) m - 2 




m/ (2 — m) 


m{m - 1] ■ {z s /\) m ~ 2 



Table C3 

Asymptotical behavior of the correlation function in velocity space. Numerical constants 

C-y m AND C-y m ARE GIVEN IN APPENDIX A. 
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The term has a finite value at a zero separation and is not important at small distances. Correlation 
at the very long distances is, however, dominated by velocity term if 7 > 2 — m. The above results are 
illustrated in Fig. |C8. 
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Fig. C8. — Correlation function in velocity space for various choice of 7 and m. Upper set of panels shown the 
density term £ p the line of sight, £ p (i2, 0), and in a transverse thin slice, £ p (Q,z s ). Asymptotics also include underlying 
r 1 slope for density correlation function before velocity mapping. 

In view of large values of A the regime r/A —>■ is the most important. 

CORRELATIONS BETWEEN DENSITY AND VELOCITY 

Our assumption in the body of the paper was that the density and velocity are uncorrelated. Therefore 
it is interesting to study whether possible velocity-density correlations alter our result. To do this we 
have to assume a statistical model for the density field. 

In what follows we shall consider the Lognormal density distribution. Particularly, we assume that 
the density can be expressed as the exponent of the Gaussian field z(x), 



p(x) = p e 



2 (x) 



(Dl) 



The advantage of this model for density is that it can describe nonlinear density fluctuations \6\ > 1 
while remaining simple^]. 
We parametrize the Gaussian distribution for z(x) in the following way 



P(z) 



1 



exp 



(2tt) 1 /M 1 /2 

34 The Gaussian model of the density field cannot describe large fluctuations and therefore we do not consider it 



z 

2A 



(D2) 
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Then the average and the second moment of the density distribution, which can be directly measured, 
are related to the parameters po and A as 



P=(p) = Poe 



A/2 



PV A 



(D3) 



In the absence of density inhomogeneities A — * 0, p — > po- 
Power spectrum in velocity space is a Fourier transform of the kernel (see eq. (0)) 

(e i/fczMl) -" z(2)] p(l)p(2)) 



(D4) 



where u z is z-component of velocity, p and 1, 2 refer to two points in space, k z is a wavevector component 
in the velocity space. 
To calculate the statistical average in eq. ( |D4|) we use the variables 



y = f/X K(l)-n 2 (2)] 
z+ = z{l) + z(2) 



(D5) 
(D6) 



which are Gaussian random quantities obeying bivariant joint Gaussian distribution. It is well known 
that a general Gaussian distribution of N correlated variables (with a zero mean values) described by 
a vector X = (xi, x 2 , ■ ■ ■ , x^) is given by the probability function 



P(x 1 ,x 2 



x N 



{2-k) n / 2 \\C\\ 1 / 2 
/ (xj) (xix 2 ) • • • (xiX N ) \ 



cxp 



-XC _1 X 



C 



(x 1 x 2 ) (xf) 

\(X!X N ) (X 2 X N ) 



• • (x 2 x N ) 
..{x%) .J 



(D7) 



We shall characterize the correlation properties of a pair (y, z + ) by three functions 

(y 2 )=D z (r) 
(4)=B(t) 
(yz+) = (u z (l)z(2)) - (u z (2)z(l)) = F(r) 



(D8) 



where F(r) is a cross correlation function^] of velocity and density, D z (r) is the z-projection of velocity 
structure function given by eq (|17D and B(r) can be expressed through the density correlation function 
£(r). Substituting those functions in eq. ( |D7| ) we get the joint distribution function 



2tt{D z B -F 2 ) 1/2 
using which one can easily compute that 

£(r) = (p(l)p(2))/p 



exp 



By 2 - 2Fyz + + D z z\ 
2(D Z B - F 2 ) 



(e z+ ) 



,B(p)/2 



(D9) 



(D10) 



and, thus, B(0) = AA = 2lna 2 /p 2 . 

35 In the equation for F(r) we have used the fact that a random vector quantity (velocity in our case) always has a 
zero correlation with any scalar quantity (as density) at the same location in space, i.e, (u z (l)p(X)) — (u z (2)p(2)) = 0. 
Indeed, the correlation functions of density-velocity in one point can be obtained by correlating velocity and density in 
two different points 1 and 2 and then bringing the points together. It is obvious that if points 1 and 2 are along z-axis 
and very close together, (u z (l)p(2)) w — (u z (2)p(l)). This proves that the correlations are zero when the points coincide. 
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Averaging kernel 



with the distribution function given by eq. ([D£]) yields 

(k z X) 2 D z (r) 



P s {K,k z ) = J d 3 re ik r exp 



exp [ik x \F(r)] £(r) 



(Dll) 



which generalizes eq. (|Tg) for a non-zero correlation of velocity and density. 

Let us now project P 3 to P 2 in a thin slice and for high |K|. Using eq. fl2"E| ) and P s given by eq. ( pil[ ) 
we obtain 



(2tt) 1 /2A 



<ir e 



KR 



Dl /2 {r) 



exp 



F 2 fr) 



2D, 



exp 



2A 2 Z^(r/A)J 



(D12) 



This expression is an analog of eq. ( p9|) in section 3 for F ^ 0. For high |K| the last exponent may be 
shown to be ~ 1 and then the equation become analogous to eq. fl4"0|). 

An upper limit for F cross-correlation function follows from general Cauchy-Swartz inequality (see 
Mathews & Walker 1970), (yz + ) 2 < (y 2 ) (z 2 + ), 



F 2 (r)< D 2 (t)B(t) = 2D 2 (t) ln£(r). 



(D13) 



and this limiting case is achieved when density (logarithm) and velocity are perfectly correlated and 
are not independent statistical quantities. If we adopt this limit, £(r) cancels out and the velocity 
fluctuations in the thin slice limit are determined by the spectrum of random velocity only: 



(2tt) 1 /2A 



(D14) 



It is easy to see that equation eq. ( pi4| ) coincides with the expression for P 2v and therefore in the case 
of of long-wave dominated density field we reproduce the thin-slice asymptotics found on the assumption 
of no correlation of velocity and density. Indeed, our analysis in section 4 showed that in this regime the 
P 2v part of the spectrum dominates the signal. Formally, we do not reproduce the result for the short- 
wave dominated density regime, where P 2p part of P 2 is important. However, we assumed the absolute 
maximum (unrealistic!) of the velocity-density correlations. Therefore, for more realistic cases we do 
expect density to reveal itself within the thin slices in the regime of short-wave density fluctuations. 

Eq. ( pilj ) also indicate that the criterion for the velocity slice to be thick may be altered if the 
velocity-density correlations are present. Apart from the criterion (|35|) one need to require that C/X > 
(1/A|K|)F. For realistically small F this condition is fulfilled when ( P5| ) is fulfilled. In the case of 
upper limit of correlation given by eq. ( pi3| ) the thickness of the slice should be increased by a factor 
of the order of a few. Asymptotics in the regime of thick slicing given by Table 1 stay the same when 
velocity-density are correlated. 

All in all, thick slice regime does not depend on velocity and density being correlated. This regime can 
provide us with the information whether density is long or short wave correlated. The density dominates 
for "very thick" slices (see criterion fl50|)) and for them it is possible to determine the density spectrum 
irrespectively of velocity- density correlations. The regime of thin slices depends on the velocity random 
field if the density field is long-wave correlated and the spectra of both velocity and density field if the 
density is short-wave correlated (see Table 1). In the case of a thin slice we recover the velocity power 
spectrum for the long-wave dominated regime. The only limiting case when velocity-density correlations 
can matter is the thin slice for the short-wave dominated density field. In the extreme case of a perfect 
velocity-density correlation the emissivity depends only on the velocity spectrum, while the emissivity 
does depend on density if the correlations are absent. This means that if the analysis of the thick sliced 
data reveals short-wave correlated density, an additional care for the analysis of the data in the thin 
slice regime might be needed. 

Our analysis of velocity-density correlations was performed for the case of the HI turbulence study 
in Galactic disc. However, it is also applicable to the studies of turbulence in individual clouds (see 
Appendix E). 
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APPLICATION TO INDIVIDUAL CLOUDS 

In the main body of the paper we were concerned with the studies of HI in the Galactic disc when the 
Galactic rotation curve served as a distance indicator. If studies of individual clouds are concerned (e.g. 
high latitude clouds) no analog of the distance-velocity relation exists. Similar problems arize when we 
deal with external galaxies. 

Here we derive the expression for 2D spectrum of HI emission applicable to studies of turbulence 
within individual clouds. To do this we modify our formalism presented in the main text to account for 
the finite size of the emitting region. We will generalize map (0): 



X S = X 



A 



f- L z - u(x) • z , (El) 



where the parameter A is just a conversion factor which specifies the units of z s coordinate. Our previous 
map (0) corresponds to the choice A = f. In the present form the map is applicable to a zero shear, 
which corresponds to f^ 1 — > 0. 

Our further treatment of the problem repeats the steps discussed in the main text but we take into 
account the finite extend of the cloud. Since in the image plane we consider only short scales relative 
to our object size, we can simplify our task by accounting for only for a finite thickness S of the cloud 
along the line-of-sight. In place of eq.(Q) the Fourier component of density in velocity space is now 

p s (k) = J* dze iAf ~ lk *- z J d 2 Xe iK -*p(x.)e- lAk * u * {x,) , (E2) 

where line-of-sight and image plane transforms are now treated separately and we use convention k = 
(K, k z ). 

To calculate the statistical average (p s (k)p*(k')) for Fourier components of the density in velocity 
space we assume statistical homogeneity for the quantities in real space, which entails that 

E 2 (k z , k' z , r) = (p(x)p(xy^ u ( x > 2 - fe >( x '> g ]) (E3) 

depends only on the difference r = x — x/. Choosing new variables r and x + = (x + x/)/2 we obtain fj|] 

(p s (kW(k')) = 25(K - K') [ S dr z f'^' 2 dz + ( d 2 RE 2 (k z , k' R, rz y^ e ^[Kz+-kir z ] (E5) 

JO Jr z /2 J 

Here fc+ = (k z + kQ/2 and k~ = k z - k' z . 

Assuming that density and velocity fields in galactic coordinates are uncorrelated (cf. Appendix D) 
and described by a correlation £(r) and structure (see eq. (]TT|)) functions we get 

E 2 (k z , k' z , r) = ^ r)e ~^[ {kz -KfD z{ oo)/2 +kz KD z{ r)] (E6) 



Important feature of eq. (|E6| ) compared to the kernel in the main text is the explicit presence of the 
velocity dispersion < u 2 z >= D z (oo)/2 in the exponent, since we have to treat k z and k' z as separate. 
This means that one cannot use pure power-law structure functions do describe velocity turbulence, 

36 Depending on convenience, we choose one of the following transformations of the double integral along a pair of 
lines-of-sight through the slice of finite thickness, z + — (z + z')/2, z~ = z — z'\ 

/■S fS rS rS-z-/2 

dz dz'F(z,z')= dz~ dz + [F{z + ,z-) + F(z+,~z-)} 

JO Jo Jo Jz-/2 



/2 

/ dz+ / dz-[F{z + ,z-) + F{S-z + ,z-)} (E4) 

Jo J-2z + 
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but the turbulence must have a maximum scale. This is perfectly in line with the introduction of finite 
cloud size along the line-of-sight. In the absence of additional physics the size S also serves as cutoff 
scale for the structure function D(r) ~ D(oo)r m /(r m + S m ), but we notice that if this law is adopted 
for Dll(t), z-component of the structure function in the case of solenoidal turbulence is 



D z (v) 
D z (oo) 



D(oo) 



D(oo) 



'■ + S m 

cs m 



1 



m/2 



i + (r/sy 



'1 — cos 



cos 6 



z/r 



(E7) 
(E8) 



For clouds introduce in Appendix F 3D spectrum in velocity space, but here use eq. (|E5| ) directly 
to obtain 2D spectrum of intensity in velocity slice from eq. (20). Expressing left-hand-side intensity 
correlation function as the Fourier integral of -P 2 (K) (see eq. (^)) and right-hand-side 3D density 
correlation in velocity space as double Fourier integral of the (p s (k)p s (k')) given by eq. (|E5| ) we can 
carry out integration over k z and k! z to obtain 



Po(K) oc / d 2 Re 



iKR 



sv 



SV-Z-/2 



dzt x 



(E9) 



[Bi(r)D z (v) 



-e 



-(/- 1 r z -2-) 2 /B(r) 



S-r z /2 
r z /2 



dz e 



+ ^-(f~ 1 z+-z+) 2 /2D z (r) 



£(R,r 2 ) 



where i?i(r) = D(oo)[l—D z (r)/2D(oo)] only weakly depends on r and ~ D(oo). Since our z s coordinates 
are in velocity units, integration over them is done directly through velocity slice of 5V width. As one 
expects, parameter A drops out of equations. Although eq. (|E9|) can be further evaluated, resulting 
error functions are not very illuminating. We shall analyze the limit f~ x — > instead, which describes 
the case when coherent motions through the cloud can be neglected. It is useful to change the order of 
integration over z+ and z~ according to (|E4fj . Then 



P 2 (K) oc / d 2 Re 



4KR 



dr z (S 



sv 

^dzt- 



7-BiW 



+ e 



-(SV-zff/B^r) 



-erf 



£(R,r 2 



(D z /2y/* 

(E10) 

The structure of this expression looks somewhat different from eq. (^3|) mainly because here we have first 
performed integration over k z , rather than z s . This form is convenient to discuss variation of the result 
with the slice thickness. We can now distinguish three regimes. If the velocity slice is thick, namely 
larger than velocity dispersion on the scale of the cloud 5V ^> B^ 2 ~ D 1 ^ 2 (oo), then an integration 
over the whole line (i.e. up to infinity) is appropriate 



Po(K) oc / d 2 Re 



iKR 



dr z (S - r z ) tan" 1 \(B l /2D z ) 1/2 } f (R, r z ) 



(Ell) 



The velocity effects disappear in this regime, with only minor residuals through the almost constant 
function tan -1 . 

When the slice is not very thick, SV <ti -D 1 ^ 2 (co), exponential terms in eq. ( |E1(J| ) are close to unity and 
one returns to the situation discussed in the paper. Namely, the slice is thin when 5V < Dl^dK.]^ 1 ). 
Then the error function in eq. ( |E10|) can be expanded in series and in the leading order 



P 2 (K)| t oc Jd 2 Re iKR J o S dr z (S-r z )B; 



1/2, 



D l J 2 (vyn 



(E12) 



This recovers the eq. (KJ) with modification for the finite size of the cloudF] . As -Bi(r) varies slowly 
with r it can be considered as a constant ~ D(oo) ~ CS m . Indeed, it is easy notice (see eq. (|E8|) ) that 

37 The presence of r z in cq. ( E12| ) results in a subdominant term that decreases faster than the main term for large |K|. 
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D z (y)/CS m = {(r/S) which plays the role of D z (r/X) in eq. fl40|). Thus the size of the cloud S plays 
the role of the correlation scale A. Also, the physically transparent condition SV < .Dy 2 (|K| _1 ) exactly 
corresponds to the transition point between thin and thick slice regimes derived earlier in flSpp . The 
fact that in a completely different regime we obtained the same condition for the separation between 
thin and thick slices shows that this condition is universal and the cancellation of / in eq. ([36|) was 
not accidental. This finding allows us to apply our theory to both individual clouds and to external 
galaxies without caring too much about a regular shear. The thick and thin slice asymptotics are given 
by Table 1 with S being used instead of A. 

The effects of the emitting region being finite is seen in the thick slice regime. This regime corresponds 
to an intermediate value of 5V, -Dy 2 (|K| -1 ) < 5V < D 1 / 2 (oo). Obviously enough, the integration over 
the full velocity dispersion provides the statistics that depends only on the density field. Inversion of 
this statistics provides the 3D density power spectrum as it is shown in L95. 

The opportunity of determining density spectrum before dealing with the velocity is very welcome. In 
other words, individual HI clouds present an excellent example for studying the interstellar statistics. 
Our results here are directly applicable to optically thin species in molecular clouds. 

ID SPECTRUM ALONG LINE OF SIGHT 

Calculation of the one dimensional power spectrum along a line-of-sight P\{k z ) in velocity space is 
very straightforward under assumptions adopted in the paper. Indeed 



or, using eq. (|Tj]) 

P x {k z ) = e - f2k * v T J dze ikzZ i(z) exp 



P 1 (k z ) = I dKP s (K,k z ) (Fl) 

(k z X) 2 D z (z/\y 



(F2) 



Asymptotics for the density P lp part (described by a correlation function £(r) oc r 7 ) and velocity term 
Pi v (which corresponds to £(r) = const) are given in the table: 



k z \ < 1 ; 


k z \ > 1 


\-\\/r f +n P lp {k z ) : (M) n+2 ; 


(£;,A) 2 ( n + 2 )/ m 


X^Pivih) ; (k z \y- m ; 


(k z \)- 2 l m 


Table F4 



Asymptotics of the ID spectrum along the line of sight at scales larger than gas sound 

SPEED, k z < l/(fv T ). 

This results are valid for 7 < 1, so n < —2. At velocity scales below gas sound speed k z > l/(fvr) 
the line-of-sight power is suppressed. Note that the expression for P\ p term in the longwave regime 
would correspond to the density spectrum along the line-of-sight in the absence of velocity mapping. 
The asymptotics are illustrated by the full solution, available for m = 1 case. 



Pi(fc 2 ) oc e^) 2 ^/^ 2 cos (n + 2)tan- 1 (l/£; z A) (^A)'(l + (^A) 2 ) v " ' ~" " (F3) 



(n+2)/2 
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To get P lv term one should substitute —3 instead of n. This spectrum allows to distinguish the short 
and long-wave density regime easily. 

For all practical purposes in the case of the Galaxy A is sufficiently large that we shall be in k z X 3> 1 
regime. In this regime the scale when the contribution from P iv and P\ p terms are equal is 

(k z \? /m = A/r (F4) 

and does not depend on n. However, which scales are dominated by the P± v term and which by the P\ p 
term does depend on whether n < —3 or n > —3 If the underlying spectrum is steep, n < —3, the short 
scales k z X > (A/ro) m ^ 2 are dominated by velocity term while the intermediate scales 1 < k z X < (A/ro) m//2 
are determined by the density. The situation is reversed for shallow n > —3 spectra. As a consequence, 
the spectrum on small scales always has shallower slope than on the intermediate ones. 

We can incorporate the effect of finite thickness of the cloud as follows. Let us model the density field 
inside cloud as p c «( x ) = P( x ) ' -^( x > S) where p(x) is statistically homogeneous random field, while F(x) 
describes the shape of the cloud. We assume the cloud is centered at the origin and has characteristic 
size S. This size, in principle, can be different in different directions, thus the vector notation. Similarly 
to eq. (H) we obtain for the power spectrum in velocity space 

(p s (kK(k)> = e-W J rf 3 re lk - r S(k, r) J d 3 x+F(x + + r/2, S)F(x + - r/2, S) . (F5) 

where x + = (x + x')/2 (compare to eq. (|J)). Note, that different Fourier components are not formally 
orthogonal when we deal with finite clouds. We did not write down the off-diagonal elements, but the 
diagonal ones, in contrast to eq. (|?p contain not a delta-finction, but a convolution of shape functions 
F. This convolution is the easiest to evaluate in Fourier space 

F 2 (r, S) = J d 3 x+F(x + + r/2, S)F(x + - r/2, S) oc J d 3 kF 2 (kS) e ikr (F6) 

ID power spectrum with finite cloud size taken into account is giving by a modified version of eq. (|F2| ) 

P x {k z ) = e- f2k ' v T J dze ikzZ £(z)exp 

The modification is insignificant for high k z > S 1 " 1 since at small scales shape function is constant and 
we can use asymptotics in Table ([Fl]). To prove this consider a spherical cloud of constant density, 
F(r/S) = const, r < 5/2; F(r/S) = 0, r > S/2. Then 

F 2 (z/S) oc l-l\z/S\ + l \z/S\ s , \z/S\<l (F8) 

This is somewhat different from a model with a finite thickness S of the cloud along the line-of-sight, 
infitie dimensions in the image plane. In this case 

F 2 {z/S) oc 1 - \z/S\, \z/S\ < 1 (F9) 

Note, that in both expressions above for z ■€ 5 1 only unity term matters and for small scale turbulence 
eq. ( [Flf) presents the 3D spectrum in velocity space. 

The ID statistics is complimentary to the 2D statistics that we dwelt upon in the main text. Applying 
ID spectral analysis to data cubes it is possible to determine whether the spectrum is long-wave or short- 
wave dominated. 



(k z X) 2 D z (z/X) 



F 2 (z,S z 



(F7) 
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